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On the Theory of Matrices. 

Bt Henry Tabee. 



I. — Elements op the Theory. 
Introductory. 

§1. Oayley, in his Memoir on the Theory of Matrices (Phil. Trans., 1858), 
defined a matrix as " a set of quantities arranged in the form of a square,"* this 
notion arising " from an abbreviated notation for a set of linear equations." 
Accordingly, Cayley laid down the laws of combination of matrices upon the 
basis of the combined effect of the matrices as operators of linear transformation 
upon a set of scalar variables or carriers. The development of the theory, as 
contained in Cayley's memoir, was the development of the consequences of these 
primary laws of combination. Before Cayley's memoir appeared, Hamilton had 
investigated the theory of such a symbol of operation as would convert three 
vectors into three linear functions of those vectors, which he called a linear 
vector operator. Such an operator is essentially identical with a matrix as 
defined by Oayley ; and some of the chief points in the theory of matrices were 
made out by Hamilton and published in his Lectures on. Quaternions (1852). They 
were, however, made out as theorems concerning linear vector operators, and 
developed by quaternion methods, through the effect of these operators upon 
vectors, and not upon the basis of the primary laws of combination above 
referred to. Nevertheless, Hamilton must be regarded as the originator of the theory 
of matrices, as he was the first to show that the symbol of a linear transformation 
might be made the subject-matter of a calculus. Cayley makes no reference 

* Oayley speaks also of rectangular matrices, and to some extent develops their theory ; he even 
alludes to the possibility of attaohipg a more general meaning to the word. His memoir deals, how- 
ever, almost exclusively with square matrices ; and as the present paper relates exclusively to such, we 
shall make no further reference to other than square matrices. 
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to Hamilton, and was of course unaware that results essentially identical with 
some of his had been obtained by Hamilton ; and, on the other hand, Hamilton's 
results related only to matrices of the third and fourth order, while Oayley's 
method was absolutely general. The identity of the two theories was explicitly 
mentioned by Clifford in a passage of his Dynamic, and was virtually recognized 
elsewhere by himself and by Tait. Sylvester carried the investigation much 
farther, developing the subject on the same basis as that which Oayley had 
adopted. Subsequent to Oayley, but previous to Sylvester, the Peirces, 
especially Charles Perrce, were led to the consideration of matrices from a dif- 
ferent point of view ; namely, from the investigation of linear associative algebras 
involving any number of linearly independent units. In this aspect, the quan- 
tities arranged in a square are looked upon as scalar coefficients of the several 
units or "vids" of an algebra belonging to a certain class. Finally, it must be 
mentioned that Hamilton and his successors made use of an interpretation of a 
linear vector operator which consists in regarding the operator or matrix as 
representing a homogeneous strain ; and in this light the axes of the strain play 
an important part in the theory. These axes may, however, be utilized without 
any reference to this interpretation, the analytical definition of the axes being 
obvious ; and I have made much use of them in the following investigation. 

This paper originated in an investigation upon the development of Clifford's 
geometrical algebras ; the consideration of the linear vector functions of these 
algebras led me to think of investigating the theory of matrices viewed as linear 
vector operators. For matrices in general these algebras furnish an instrument 
analogous to that which quaternions affords for the investigation of mat- 
rices of the third order. I shall, however, reserve to a subsequent paper the 
consideration of this particular system of algebras (of which I have obtained a 
tolerably complete development) and its utilization for the theory of matrices. 
In the present paper I regard a matrix as a linear unit operator, operating upon 
the linearly independent units of an algebra, without reference to any meaning 
of such units, or to any properties which these units may have in combination 
with each other ; and I have in this way endeavored to develop the theory of 
matrices. From this point of view I am able to give a very simple proof of 
Cayley's "identical equation," and also of Sylvester's most important results, 
the law of latency, the law of nullity, and Sylvester's formula for any function of 
a single matrix. I have also completed the investigation of the nullity of the 
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factors of the identical equation (Sylvester's "corollary of the law of nullity") 
which Sylvester had treated only in the special case when all the latent roots of 
the matrix were distinct ; and I have shown that in addition to nonions there are 
an infinity of algebras exactly analogous to quaternions. This analogy I have 
also extended, and it appears that every matrix is resolvable, precisely as a qua- 
ternion is, into the product of a tensor and a versor; the latter gives rise in a 
matrix of any order to functions analogous to the sine and cosine to which 
these functions reduce when a — 2 ; and thus I find that there exists an infinity 
of formulae analogous to De Moivre's. Finally, I have shown that the whole 
theory of matrices may be regarded as contained in the theory of Clifford's 
geometric algebras, i. e. in the theory of sets of quaternion algebras which are 
such that the units of one set are commutative with those of any other. Other 
results I have obtained not immediately connected with the pure theory of 
matrices, but having reference to the matrix viewed as a linear unit operator. 

These results are contained in the second part of this paper ; the first part 
contains only the elementary notions and theorems developed from the point of 
view of the matrix as an operator, and is not necessary to the understanding of 
the second part, at all events by one acquainted with Cayley's memoir. To this 
must be excepted §9, which contains an account of Charles Peirce's system of 
quadrate algebras and their connection with the theory of matrices. 

I have in. I §10 given a slight sketch of the history of the theory. 

Definition of a Linear Unit Operator. 

§2. The extension of any selection of a of the units of any algebra will be 
termed a ground of order w. A linear unit operator q> of a ground of order a is 
an operator which converts each of a linearly independent quantities in the 
ground into a quantity in the same ground, and which is such that 

<f> (ma + nn) — m ((pa) + n (tyn) , 

where G and n are any two quantities in the ground, and m and n are scalars.* 
Thus, if a lt a % , . . . . a u are a linearly independent quantities of a ground of 

order a, and if 

p = x x a x + £Ca« 2 +....+ » w a w . 

* Of course this requirement is equivalent to the requirement that ^ shall be distributive over any 
sum of quantities in the ground ; the equation would then follow at once for m and n positive integers, 
and thence very simply for m and n any scalars. 

44 
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then $p = x x . ^>ft x + x% • 4>«2 +••••+ s»o> •4 )a » > 

where, since the 4>a's are in the same ground as the a's, and these are o in 
number and linearly independent, we may put 

# a i = «n«i + «2i«2 + ••••+ ««,ia w , 
$a 2 = a 12 a x + a 22 a 2 +....+ a u2 a M , 



whence the linear unit function typ becomes equal to 

(a u X! + a n Xi + + a loi x a ) a x + (a^ + a 22 z 2 + + a 2w a; w ) a a + . . . . 

+ («„ia;i + a a &% + ....+ a aa xj) a w . 

By choosing the a's as the units of the set (x lt x z , . . . . x a ), it may be 
regarded as identical with p, when it is evident that <£> is the matrix 

( an a iz • • • • <*!„ ) 
^21 ^22 • • • • a% m 



Thus, either the successive columns of the array constituting the matrix may be 
regarded as operating upon the successive quantities (a x , a 2 , . . . . a„) represent- 
ing the ground, or the successive rows may be regarded as operating upon the 

coefficient's (a^, x % a;„) of these quantities. 

Regarding <p as operating upon the set p = (a^., £C 2 , . . . . #„), i. e. upon the 
coefficients of the a's in the expression for- p , the resulting set consisting of the 
coefficients in the expression for the linear unit function <£p may be denoted by 

((<HO> (4>|ai 2 ), (<M »„)); and so 

$p = (4>|x 1 )a 1 + (<j!>|a: 2 )a 2 + + ($ | x„) a„ . 

In virtue of the identification of a matrix and a linear unit operator, these terms 
will hereafter be used interchangeably. 

The coefficients of the a's in the expression for the <£>a's which form the 
columns of the matrix, will be termed its constituents. 

The converse or transverse of <p, denoted by <j!>, will be defined as the matrix 
formed from $ by interchanging the rows and columns of its constituents.* Whence 

• The term transverse was used by Cayley in his memoir to denote this function of a matrix, and 
has subsequently been used by Sylvester in the same sense ; the term converse was that employed by 
Charles Peirce for this purpose, and the symbol of conversion or transversion employed here is also his. 
Hamilton used the term conjugate of $ to denote the converse of <j> , and denoted it by <!>'. 
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where 



$a 2 = a u a! + a 12 a a + ....+ a 1<0 a w , 
$a 2 = a 21 ax + <*22«2 + ....+ a ia a a , 
etc. 



The relation " converse of" is evidently reciprocal, thus ^ = <p. 
Two matrices of the same ground are equal if they have invariably the same 
effect upon the same quantity. Thus if 4> is the matrix 

( hi ha ■ • • • ha, ) 
hi ht • • • • hi 



hi h2 • • • • % 

and the a's still represent the ground, so that 

■4*ti = hi<*i + hi<H + • • 
4<a 2 — JiaOCi + 5 22 a 2 + . . 
etc., 



• + hia», 

• T" hi a ia I 



and if ^prr^p, for any quantity in the ground, then a u = 6 11 , a 12 = 6 12 , etc. 
Consequently, if two matrices are equal, every constituent of the one is equal to 
the corresponding constituent of the other ; and so an equation between matrices 
gives rise to a 2 scalar equations. 

In defining the linear unit operator <£> by its effect upon a linearly indepen- 
dent quantities in a certain ground of Order a, the law of multiplication of the 
units of the ground is not in any way involved. Whence it follows that there is 
only an apparent loss of generality in taking the units of the ground from any 
particular algebra. If the elementary units of Clifford's o-way algebra* are 
chosen, any quantity in the ground is a vector in w-dimensional space ; <£> may 
then be termed a linear vector operator, and has a definite geometrical signification, 
namely, it represents a homogeneous strain. For if p and a are vectors to two 
parallel straight lines, then we may put p = a 1 + x/3, <y = a z + y@, and by 
definition <£>p = q>a x + x<p@ , 4><r = 4><x 2 + ytyfi, so the displacement of points in 
space which 4> effects is such that parallel lines remain parallel ; whence, exten- 
sions of any order which are parallel remain parallel after the application of the 



* Clifford's w'way algebra is an algebra linear in the product of even order of o> " elementary units," 
ilia .... i» , such that ii> = — a, and i? = — 1 . See this Journal, Vol. I. 
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strain <p. Without, however, determining the nature of the units of the ground, 
i. e. the properties its units may have in combination, two extensions belonging 
to the ground of order m and n respectively, where rn<Cn, may be termed 
parallel if every quantity in the one may be represented, for some set of values 
of the x's, by 

and every quantity of the other, for some set of values of the y's, by 

<*2 + yift + 1/A + + 2/m/3m + + VnPn, 

the a's and /3's being quantities of the ground. A homogeneous strain may then 
be defined as a displacement in a certain extension (namely, that of the ground) 
which leaves unaltered the parallelism of any two extensions of the ground. 
Obviously 4> effects such a displacement of quantities. 

A consequence of some interest follows from the identification of the theory 
of linear vector operators with the theory of matrices : as quaternions is identical 
with the theory of dual matrices, and thus with the theory of homogeneous 
strains in a plane, to every proposition concerning space of three dimensions (or 
of four dimensions) which can be proved by quaternions, corresponds a proposi- 
tion concerning the kinematics of a plane, such that when either is proved so also 
is the other. 

Hereafter it should be understood, that when quantities are spoken of as 
the operands of a matrix <p, they are to be regarded as in the ground pertaining 
to 4> even when this is not explicitly stated. 

Addition of Matrices, Multiplication by a Scalar, and Form of a Scalar. 

§3. Employing the notation of the last section for q> and 4', let ^ denote 
the matrix 

( <h\ + hi «ia + &12 «i» + ^lo, ) 

a 21 "T" ^21 a 22 H" 0»S • • • • a 2a + Ofy 

etc., etc. 
It is obvious from the definition of a linear unit operator that for any quantity p , 

Z? = $? + ty- 
Whence this sum may be denoted by (<$> + 4*) p, giving the equation in matrices, 
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( «u 
etc. 



a u 



•«!•) + ( &n 



«2a> 



hi 

etc. 



?>22 



&la, ) = ( «11 + &11 «12 + &tt «ta + &ta ) 

• • C*2a, H~ 2a) 



«21 + hi «2S + &2: 

etc. 



As a consequence of the nature of scalar addition, it follows that the addition of 
matrices is commutative and associative, and such that its inverse, subtraction, is 
determinative, or may be regarded as the direct operation performed with an 
inverse quantity, the negative. The negative of <£>, which may be denoted by 
( — <£>), is of course that matrix whose constituents are the negative of the 
corresponding constituents of <J!>; otherwise <?>p+( — <?>)p would not be zero for 
all values of p. 

Evidently the converse of the sum of two matrices is the sum of their con- 
verses. Thus : . -~r^ , . -,~t 

§4. From the last section it follows that if to is any positive integer, 
(rn\ an a 12 . . . . a la ) = ( ma n ma 12 .... rna loi ) 



a 21 a 22 

etc. 



• a Zo> 



ma n ma z2 
etc. 



ma 



2«> 



Thence it may be shown very simply that the equation holds for any scalar 
to. And hence, obviously, for any scalar to , (m<p) p = q> (w?p) = to (<pp) . If, then, 
n is any other scalar, for any quantity p in the ground we have (m.nq>)p 
= m(nq>.p) = to (n.typ) = rnn.typ = (mn.<p)p; also (to + n)$.p = (m + n).typ 
= to .typ + n . typ = to<£> . p + nty . p — {rnty + nty) p ; consequently, m.nty = mn .<f> 
and (to + re) <£> = to^> + to£, and therefore scalars as multipliers of matrices have 
all the properties they possess in combination with monomial or scalar quanti- 
ties. Hence in future the combination rn.(pp = vnty.p may be denoted by m$p, 
etc. 

A scalar is itself obviously a linear unit operator or matrix ; but the only 
matrix which merely multiplies a set or quantity by a scalar to is 

(toO etc. ) 

to etc. 

to etc. 
etc. 



where all the constituents are zero except those in the principal diagonal, and 
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they are all equal torn. If m = 1 , this array gives the matrix form of unity ; 
if m = 0, it gives the matrix form of zero. It will presently be shown that the 
product of two matrices is also a matrix, and that the product is distributive 
over the operand. In conjunction with the matrix form of a scalar, this makes 
significant the compound 4>m , m, being a scalar. For since for any quantity p 
(<pm) p = <£> (mp) = m (<fy>) ; hence q>m = mty . 

Multiplication of Matrices and its Inverse. 

§5. The combination 4>i|f must be defined as the result of operating by the 
linear unit operator q> upon the linear vector function or quantity $p, i. e. ^'P 
= ^> (ij'p) . With the same notation as before, for $ and $, let 

% = ( Z r a lr b rl A r ai r o r2 .... z r a lr b r(a ) 



2* r ci wr o r i 2^ r a ur b r 2 .... 2< r a ar b r 



Then for any quantity p , 

<Nf = 4> PA a • «i + 2 A A • «a + • • 

-— \ a u • Zgbi s x s -f- a^.^gb^gXg -J- . . 

+ («21 • 2 A A + «22 • 2 A A + • • 
— Zi r 2j s (l] r rs X s . Otj -f- ^/ r ^/ a a^ r O ra X a . CLo. -p . . 

= (2 r a lr b rl . Xi + 2 r a lr b ri . as 2 + . . 

+ (2 r «2 Al • «1 + 2r«2 A2 . #2 + • • 

Whence we may put ^ =t ^ , or 

( on «i 2 . . . . a lia J b n b 12 . . . . b u ) 



• + ^gb^Xg . <x M J 

. + a la .2gb 01 gX,)a 1 

. + a^.Xgb^Xg) a 2 + etc. 

. ~J~ 2t r 2/ s (l m -0 rs X s . GC.„ 

. + 2 r a 2r &, w .a; w ) a 2 + etc. 



a 21 <% 



' ^2<a 



bo, b 



'21 "22 



^0)1^02 • • • • <^» w 6<ol6, 

( «n in + «i2 &21 + • • • 

a 21 &11 + a 22 2 i + . . . 



• + «i» &a,D «n &12 + «i2 ^22 + • • • • + a i». Ki, etc. ) 

. -j- a 2o> "<ol j a 21 #12 "T ^22 "22 "T • • • • "T" #2a) "w2 1 ©tC. 



a <ol"ll *T ^2 "21 ~T • • • • ~T" a w»"a)l) <*«)lOl2 + a <o2"22 +•••• + a »J>tata ! e tC 

i. e. the product of the two matrices 4) and ^ in the order named is formed as 
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follows : The constituent of the r th row and s th column of <pTp is the sum of the 

products obtained by multiplying each constituent of the r th row of $ into the 

corresponding constituent of the s th column of 4'- 

Since by definition the product of two matrices is associative with their 

operand, it follows that the multiplication of matrices is also associative. For 

by definition 

<?> fix) ■? — ♦• ("h) P = 4> W-zp)> 
and likewise, 

(tf>) X -9 ~ (^) •ZP — VW'-Z?)- 
As this is true for any quantity in the ground, hence 

Q&x) = (#)%• 

The multiplication of matrices is also distributive over addition. For if the 
quantities o' and o" denote severally the linear unit functions <pp and 4p of any 
quantity o in the ground, by §2 and §3, 

4> (4 + X) P = <?> Wf + %p) = * (p' + p") 
= 4>p' + 4>p" = 'H'P + 4>%p 
= (^4 + 4>%)p- 

The commutative principle does not in general hold in the calculus of matrices. 

The converse of the product of two matrices is the product of their con- 
verses in the reverse order, as may readily be proved, This gives the formula 

§6. Division is the operation inverse to multiplication. Since multiplication 
is not in general commutative, two signs are required for division. In the last 
section, two matrices being given, it was required to find their product in either 
order ; the problem inverse to this is to find a matrix which when multiplied 
into or by a given matrix <£> shall have as product a given matrix -<l>. In general, 
this problem is susceptible of one or more solutions : the matrix or matrices 
which when multiplied into <£> give ^ as product may be denoted by (^ :<£>); and 

those which when multiplied by $ give 4> as product may be denoted by — or -IT . 

Division is therefore defined by the equations 

(4-:<?>) $ = 4>, 



♦(£) = +• 
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Invertible multiplication is multiplication whose inverse is determinative. 
If the multiplication of matrices is invertible, it would follow that if %$ = 4'. 
and ^ = '4', then £ = #'; likewise, if <p% = 4' and <£>£' = 4*, then % = %'. 
Hence division would be subject further to the conditions 

by means of which the above results may be obtained immediately. To assume 
that multiplication is invertible is equivalent to regarding the inverse operation 
as the direct operation performed with an inverse quantity or matrix, which 
shares the associative and other properties of matrices in general. This inverse 
matrix is termed the reciprocal.* 

It may now be shown that in general there exists a reciprocal, and hence 
that, in general, multiplication is invertible. From the rule for the formation 
of the product of two matrices, it follows that, if from the array of constituents 
representing a matrix q> another matrix <3> be formed, in which each constituent 
of the first array is replaced by the logarithmic differential derivative with 
respect to that constituent of the determinant of the array (provided this deter- 
minant is not null), the product of q> and the converse of <£ in either order is 
equal to unity, f Hence <E> may be denoted by qr 1 and is termed the reciprocal 
of #. In other words, if |^| denote the determinant of the array representing 
the matrix (which will in future be termed the content of the matrix), and if 
A^ represent the differential operator 

( _3_ _d_ _d_ ) 

da n da 1% ' da u 

d d d 



da™ da«, 



da, 



2u> 




* The term invertible multiplication was employed by Mr. C. S. Peirce in his Logic of Relatives, 
Memoirs Am. Acad., Vol. IX, in which, and in Mr. Peirce 's other writings, will be found the substance 
of this account of division. Mr. Peirce has shown that the only algebras in which division is always 
determinative are ordinary algebra, with and without the imaginary, and real (semi) quaternions, this 
Journal, Vol. IV. 

t In the wording of this statement of the relation between $ and * I have followed Prof. Sylvester, 
this Journal, Vol. VI. 
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then ^_j _ A^\ 



1*1 

A matrix whose content is zero is termed vacuous. The inverse processes 
are sometimes possible with such a matrix ; but since a vacuous matrix has no 
reciprocal, the results are indeterminate. Thus if 

m— ( a ) 

| ma I ' 

% = {b 0) 

\ mb |'' 

the matrix (-%- J which when multiplied by m gives £ as product, is 

(A > 

a 

c d 

where c and d are any two scalars. 

If $£ = 0, it is evident that either both q> and ^ are vacuous, or one or 
both are zero. 

From the definition of the reciprocal it follows that (<j5) _1 = (qr 1 ). For 

taking the converse of both sides of the equation tyqr 1 = 1 , we get (^> -1 ) $ = 1; 
hence (qr 1 ) = ($) _1 . 

§7. Regarding a matrix as an operator, the problem inverse to that of 
finding the effect of a matrix <£ upon any quantity p in the ground (i. e. of find- 
ing the product of q> into p) is, given a quantity <y in the ground, to find another 
quantity which <p will transform into <r. Only one sign is needed for the inverse 
of functional multiplication, which is defined by the equation 

(7 

The quotient — is always determinate if q> is non-vacuous. If the matrix 

<£> transforms any extension (a a , a 3 , . . . a m ) into one of lower order (0 lt /? 2 , . . . (3 n ) , 
and if a = ^/^ + ?/ 2 /? 2 +.... + y n (3 n , the problem to find a quantity which <£> 

transforms into ff is in each case possible, but the quotient — is not in each case 

45 
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determinate. If p = x&x + ayx 2 + ....+ x m a m , there does not exist in every 
case a quotient -*- . 

Powers and Roots of a Matrix. 
§8. If m is an integer, $ m is of course that matrix which results from mul- 

m 

tiplying $ m times by itself. If n is also an integer, 4>« must be defined as that 
matrix whose n tu power is the <m th power of <£ . The irrational scalar power <p p 
of <?> must be regarded, as in common algebra, as the limit of the series <p p ', <p p ", 
etc., where p', p", etc., are successive approximations to^>, i. e. if 

$p' = ( a' n a{ 2 etc. ) 

Otjji ^22 GuC. 

etc. 

and if the limits of the a n s as p f approaches the limits are respectively the simi- 
larly situated constituents of the matrix 

4- = ( «n <h» etc. ) 

^21 ^22 G vC» 

etc. 

then <p p = $. Whence q> p — ^p' may be made as nearly equal to zero as we 
please by taking p' sufficiently near to p ; and consequently, if <p p ' is susceptible 
of having a reciprocal, <p p :<p p> may be made as near unity as we please by 
taking p' sufficiently near to p. It is obvious that if <y is the limit for any 
quantity p of ty p 'p, as p' approaches the limit p, then <£ p p = <?. 

In general, the matrix $ m where m is an integer has a reciprocal which 
may be denoted by (4> m ) -1 . Prof. Sylvester shows, virtually as follows, that 
(4> m ) -1 = ($ -1 ) ro . Since <?>4> -1 = 4> _1 4> = 1 , hence <p and qr l are commutative ; con- 
sequently 1 = (4T Vj>)" = for 1 )" 1 *", or (4." 8 )- 1 = fo- 1 )^" (cp™)- 1 = far 1 )". Frac- 
tional and irrational powers of a matrix also have in general a reciprocal ; and 
in this case also it may be shown that (4> m ) J_1 = <p~ m . Whence it follows that 
for any two scalars ($")" = $ {mn \ 

It is obvious that ty m q> n = q> m+n for any two integers m and n, being merely 
an expression of the associative principle ; and it is sufficiently plain that with 
a proper understanding the equation holds for any two scalars. 
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Linear-form Representation of a Matrix. 
§9. Let the set of forms or vids* 



(o,:oi) (a, 
(a 8 :ax) (a g 



ocg) (a x : a 3 ) etc., 
a 2 ) (a 2 :a 3 ) etc., 
ag) (a 3 :a 3 ) etc., 



as operators upon the quantities a x , % , etc., be defined by the equations 

(a r :a,)a s = a r , 
(a r :a s )a t = 0; 

and, when members of the set are operands, by 

(a r : a,) (a, : a„) = (a r : a n ) , 
(a r :a„) (a«:a n ) = 0. 

Then the expression 

<J> = « u (a x : a x ) + a 12 (a x : Og) + a 13 (a a : a 3 ) + etc. 

+ <h\ ( a 3 : «i) + a n («z : a 2) + a 23 ( a 2 : a 3) + etc « 

+ a 31 (a s :a!) + a 32 (a 3 :a 2 ) + a 33 (a 3 :a 3 ) + etc. + etc., 

linear in these vids, considered as an operator upon any quantity 

p = x x a x + sc 2 a 2 + x 3 a 3 + etc., 

is identical with the matrix 

q> = ( a n a 12 a 13 etc. ) 
ct 2 i ot 22 ot 23 etc. 
a 3 i a 32 a 33 etc. 
etc. 

Moreover, if 'P is a similar expression in the vids (<x r :a s ) in which their coeffi- 
cients the a's have been replaced by b's with corresponding subscripts, then 
m$ + nf , where m and n are scalars, regarded as an operator upon p , is clearly 
identical with the matrix rnty-\-n$. Likewise <J>*P has the same effect upon 
any operand p as the matrix 4>i^; for 

* The term vid was introduced by Mr. C. S. Peirce to denote the units or letters of an algebra. It 
will be employed in what follows to denote the forms (a, : a,) , which will also be termed the elementary 
units of a matrix. 
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which corresponds to the matrix 

{ 2t r (%i r o r i 2i r (ii r o r % 2t r cti r o r g etc. j 

^r^ar^Vl •£t r Q'2rOr% ^r^ir^r3 6tC, 

etc. 

Any complete system of o 2 of these vids forms a pure algebra of a certain 
class termed by Clifford quadrates; and expressed in terms of these units is in 
what may be termed its canonical form.* I shall therefore call an algebra linear 
in g> 2 of these vids a quadrate algebra of order a ; and any expression linear in 
the vids, a quadrate form. The multiplication tables to which these algebras 
give rise are similar, and are immediately obtained from the laws to which the 
vids are subject. Thus if a= 2, let 



h 



denote a complete set of four of these vids. These letters or units give an 
algebra whose multiplication table is 



h i 



i 


3 










k 





i 


3 


I 














& 


I 



This is the algebra (<? 4 ) of Prof. Peirce's linear associative algebras, and is a form 

of quaternions. If a = 3 , let 

% j h 

I rn n 

p q r 

denote the complete set of vids of the quadrate algebra of order three ; these 

* It seems appropriate to have a term to express' that form of an algebra in which its units are 
capable of a classification, according to the requirements of the analysis of Peirce's Linear Associative 
Algebra. For this purpose I employ the term canonical form. Thus the form of quaternions given 
below is its canonical form ; Hamilton's units are expressions linear in those given below. 



Taber : On the Theory of Matrices. 



351 



letters give the algebra which was termed nonions by Clifford ; its multiplication 
table is 

ijklmnpqr 
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In virtue of the correspondence that has been shown in the first part of this 
section to exist between quadrate forms and matrices, it follows that to any 
function of one or more matrices corresponds a quadrate form which is the same 
function of one or more quadrate forms corresponding to the several matrices, 
and which has upon any operand p the same effect as the resulting matrix. Hence 
whatever equality subsists between combinations of matrices, if it can receive 
interpretation as like operations upon the ground, also subsists between the 
same combinations of corresponding quadrate forms. Thus it appears that there 
is no essential difference between the theory of matrices and the theory of quad- 
rate forms. Viewed in this aspect, the scalar quantities arranged in a square 
and forming a matrix may be regarded as the scalar coefficients of the several 
vids. Or the substantial identity between the theory of matrices and of quad- 
rate algebras may be brought out by considering all the possible a 2 matrices of 
order a which can be formed with one constituent unity and the remainder zero ; 
when it is very readily seen that these matrices have the same multiplication 
table, and the same effect upon the ground as the a? vids of the quadrate algebra 
of order a ; and since any matrix of order a may be regarded as an expression 
in an algebra whose Units are these a 2 matrices, we thus have two algebras whose 
multiplication table is the same, and which, consequently, are identical. Hence 
a matrix is a quadrate form, and conversely. It is easily shown that the matrix 
of order a whose non-zero constituent is in the r th row and s th column is identical 
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with the vid occupying the same place in the quadrate system of order a. Thus 
if a = 2, the four matrices 

(1 0), (0 1), (0 0), (0 0), 

loot I 1 1 1 ol loil 



in terms of which any dual matrix may be expressed, are respectively identical 
with the units i,j, h, I of the quadrate algebra of order four, or their equals the 
vids (%:%), (a^a-j), (a 2 :ai), (a a :a.j); the matrices have, moreover, the same effect 
upon the ground (a x , a 2 ) as these vids, etc. 

The discovery of these systems of vids or of quadrate algebras is due to Mr. 
C. S. Peirce, by whom, through' this discovery, the relation of the theory of 
quadrate forms (or the theory of matrices) to the general theory of multiple 
algebra was first made clear. 

By § 3, $ = a u («! : a t ) + o» (% : a 2 ) + etc. 

+ #21 («2 : a i) + a 22 («2 : <h) + etc - + etc. 

But by definition cj> is obtained from q> by interchanging its row and columns. 
Hence 

(a r : a r ) = (a r :a r ) , (a r : a s ) — (a 3 : a r ) . 

The vids of the type (a r :a r ) I shall term self- transverse or self-converse vids ; 
they may also be termed symmetric vids. Those of the type (a r :a s ) I shall 
term non-symmetric vids. Mr. Peirce terms the vids of the first type self-vids ; 
the vids of the second type he terms alio-vids. 

Sketch of the History of the Development of the Theory of Matrices. 

1 10. An outline of the origin of the theory of matrices was given in §1. 
Subsequent to Oayley's memoir, the next advance was made in 1870 by Charles S. 
Peirce, who, in his investigations upon the extension of the Boolian calculus to 
the logic of relatives,* came upon a set of forms (considered in §9) constituting a 
system virtually identical with the calculus of matrices. Peirce showed that any 
relative term involving not more than one correlate (dual relative)f could be 
represented as an expression linear in the units of a linear transformation. 

* Description of a Notation for the Logic of Relatives, Memoirs Am. Acad. Sciences, Vol. IX (1870). 
t Such as " lover of," " loved by," " mother of," etc. ; but not " buyer of— from— ," etc. 
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Whence follows the remarkable theorem that every such relation between any 
group of objects can be represented by a matrix. As has been stated, the relation 
of the theory of matrices, as. algebras of a certain class (see §9), to linear associ- 
ative algebra in general, was first made clear through the light thrown on the 
subject by Peirce's systems of vids. 

Charles Peirce has made the great discovery that the whole theory of linear 
associative algebra is included in the theory of matrices. He has shown that 
every linear associative algebra has a relative form, i. e. its units may be 
expressed linearly in terms of the vids (denoted in his notation by (A: A), (A:B), 
etc.) of a linear transformation ; and consequently, that any expression in the 
algebra can be represented by a matrix.* Whence the theory of all possible 
linear associative algebras is only the theory of all possible sets of matrices con- 
stituting a group in Benjamin Peirce's sense, i.e. which are such that the product 
of any two members of the set can be expressed linearly in terms of itself and the 
other members of the set alone. Charles Peirce has, moreover, given the relative 
or matrix form of all the algebras considered by his father in his Linear Associ- 
ative Algebra. 

To Charles Peirce, in conjunction with his father, the identification of 
quaternions with the quadrate algebra of order two (i. e. the algebra of dual 
matrices) is also due. Cayley, in his memoir, had remarked upon the simi- 
larity a certain system of three dual matrices had to the i,j, h of quaternion ; 
but the identification was not completed until the remarkable discovery by Ben- 
jamin Peirce of a form of quaternions, which, in §9, I have termed the canonical 
form, and which results from choosing the linear functions 



1+W — 1 j + W — 1 — j+W — 1 1 — iV — l 



2 



' o ' 



of unity and any three mutually normal unit vectors i,j,k, as units of the 
algebra. These units obviously have the same multiplication table as the vids 
of a dual matrix. In his memoir of 1870 Charles Peirce had given, as an example 
of the infinite system of quadrate algebras, the multiplication table of the 
quadrate algebra next in order after quaternions, afterwards named nonions 
by Clifford ; and he states in the Johns Hopkins Univ. Circs. No. 22 (April, 

*See Mr. Peirce's Logic of Relatives, above referred to ; also this Journal, Vol. IV, p. 221, and the 
Prog. Am. Acad, of Arts and Sciences for 1875. 
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188 2), that the identification of quaternions with the quadrate algebra of order 
two suggested to his father and himself that, as nonions was thus shown to be 
the exact analogue of quaternions, there ought to be a form of nonions analo- 
gous to Hamilton's tetranomial form of quaternions, which form of nonions, Mr. 
Peirce states, either his father or himself found.* By means of this form any 
expression in the algebra is susceptible of representation linearly in terms of 
unity and eight non-scalar cube roots of unity, just as any expression in quater- 
nions is susceptible of representation in terms of unity, and three square roots of 
unity (or of — 1) . This form of nonions I shall term the octanomial form. 
The Peirces' discovery of the octanomial form of nonions was not published. 
The priority of publication of this form belongs to Sylvester, who discovered 
it subsequently to the Peirces without any knowledge of their investigations 
upon nonions. f 

To Sylvester we owe most of the development of the theory of matrices. 
In his unfinished memoir in this Journal,^ Sylvester distinguished between the 
different degree of vacuity and nullity of a matrix, replacing by these terms the 
term indeterminate used rather vaguely by Oayley to denote a null or vacuous 
matrix. It should be stated that Clifford had previously distinguished between 
the different degrees of nullity, employing the term indeterminate, with the prefix 
singly, doubly, etc. In this memoir Sylvester showed also how to derive the chain 
of equations from the identical equation, and the relation of these to the latent 
function of two or more matrices taken in a certain order. In a series of papers 
in the Johns Hopkins Univ. Circulars and in the Phil. Mag., Sylvester added 
largely to the theory of matrices.|| He demonstrated the extension of Hamilton's 
theorem concerning the cubic equation to which every matrix of the third order 
is subject, which was enunciated without proof by Cayley in 1858 ; among the 
other important results are what I term the second branch of the law of nullity, 
the corollary of the law of nullity, the law of latency, and the expression for 
any function of a single matrix. Sylvester has also extended the solution of 

* " So much was published by me in 1870 [the multiplication table of the canonical form of nonions, 
etc.], and it then occurred to my father or to me (probably in conversing together) that since this 
algebra was thus shown (through his form of quaternions) to be the strict analogue of quaternions, 
there ought to exist a form of it analogous to Hamilton's standard tetranomial form. That form either 
he or I certainly found. I cannot remember which after so many years ; whichever did must have 
found it at once." — Johns Hopkins Univ. Circs. No. 33 (April, 1883). 

t Johns Hopkins Univ. Circs. Nos. 15 and 17 (1883). 

J " Lectures on Universal Multiple Algebra," this Journal, Vol. VI (1883). 

|| Nos. 15, 17, 37, 38, 33. 
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equations in quaternions or dual matrices beyond the point where Hamilton left 
the subject. 

Clifford's contributions to the theory of matrices have been left to the 
last, though in point of time they preceded Sylvester's researches. They are 
contained in his "Fragment on Matrices," published in the posthumous edition 
of his works. I shall show subsequently that the main proposition of which he 
gives a proof is false. But the basis of his treatment of the subject is an important 
contribution to the theory of matrices. It is the method which .is adopted in 
this paper ; and the demonstration which I have given in the second part of this 
paper of the law of nullity is based on a result contained in Clifford's "Frag- 
ment." 

II. — Development of the Theory by Means of the Axes of a Matrix. 

Axes of a Matrix, or Quantities for which (<£> — g) p = , where g is a Scalar. 
§1. If p = Xx-jOii and (<£> — g) p = , then 

(<% — g) x i + a n x z + ••••+ ai«.a>a, = , 
a n x i + («22 — g) x % + • • • • ■+ a 2A = , 
etc., etc. 

The resultant of this system is called the latent function of <£> , and will be denoted 
by \q> — g\, as it is the determinant of the matrix <|> — g. From its vanishing 
results an equation of order a in g, 

g°' — in a _ 1 g'-- 1 + m m _ 2 g'— 2 — T^± m=0. 

The constant term m is evidently the determinant of the array of constituents 
forming the matrix, and has been denoted by | <p \ ; m 1 is the sum of all the prin- 
cipal first minors of |$|;' and, generally, m a _ K is the sum of all the principal x tu 
minors of |#|. The term latent function is due to Prof. Sylvester, and the 

roots gi, g z g» of the a ie he terms the latent roots of 4> . If | $ | = , it has 

already been stated that <?> is then termed vacuous, and evidently has one latent 
root zero. If, in addition, nt x ^<), q> has only one latent root zero, and is then 
said to have the vacuity one, or is simply vacuous. More generally, if all the 
w?'s from m to m K _ 1 are zero, and rn K =f= 0, the matrix <£ has just x latent roots 
zero, and is said to have the vacuity %. If \$\ =f= 0, $ is non-vacuous, or its 
vacuity is zero. 
46 
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Corresponding respectively to each latent root of 4> are the a quantities 
Pii fa, • • - • f>») suc h that 

— (<?> — 9i)pi={$ — 9%)?%— ■ ■ • •=(* — 0,)P-- 

These quantitities will be termed the axes of ' q>. 

§2. There is no linear relation between any set of the axes of <^> correspond- 
ing to distinct latent roots. For suppose any set of n axes are linearly related 
which correspond to the distinct latent roots g lt g it . . . . g n . Let % be the 
number of these axes which are linearly independent, and suppose they corre- 
spond to the first x latent roots. Then any axis corresponding to any of the 
remaining n — x latent roots can be expressed linearly in terms of these. Thus 

P*+i = fyi + t*Pz + • • • • + t$ K , 
•'• 9 , « + iPk+i = <?>P*+i= hg\?i + hg%?% + • • • -+t K g K p K ; 

but, since the tfs are not all zero, this is impossible. 

Whence, if all the latent roots are distinct, the axes of <p are all linearly 
independent. If a set of latent roots become equal, linear relations may arise 
between the set of axes corresponding to them, i. e. certain of these axes may be 
projected into the extension of the remaining axes corresponding to that set of 
latent roots, or all the axes of the set may become coincident. 

If two or more axes of the set remain linearly independent when the set of 
latent roots become equal, these axes and also the remaining axes become inde- 
terminate. Thus if the n latent roots g lt g 2 , . . . . g n ultimately become equal, 
of the axes corresponding to them only the first x may remain linearly inde- 
pendent, and the remaining n — x will then be expressible linearly in terms of 
them. These x axes p lt p 2 , . . . . p K will all satisfy the equation ($ — <ft)p = 0, 
and consequently any expression ccjp! •+■ x % p % + ....+ x K p K linear in them will also 
satisfy this equation, and hence will be an axis of <?> . In this case any x quantities 
giving the extension of p u p 2 , . . . . p K , together with any n — x other, quantities in 
their extension, may be regarded as the axes of <p corresponding to the %-fold 
latent root g v No quantity in the extension of p x , p z , . . . . p K can be in any 
linear relation with axes corresponding to latent roots other than g x by what has 
just been proved, as such quantities are axes of $ corresponding to g x . 

The matrix ^ = 4> — <7i evidently has as latent roots g x — g x , g 2 — gi, g s — <fr , 
etc.; for if g is any latent root of <p , then | 4> — (g — g t ) \ = \ (<p — g x ) — (g — g x ) \ 
= \(f) — g\ = ; and hence g — g x is a latent root of 4^ since it is a root of the 
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latent function of <^. Conversely, if g — g x is a latent root of $, g is a latent root 
of <?>. Now if the axis of q> corresponding to g x (i. e, the axes of 4> corresponding 
to zero) is indeterminate, then every first minor of \<p — grj =: |ij>| is zero; 
consequently every principal first minor of 1 4> | is zero ; but then $ has two 
latent roots zero, and consequently <?> has two latent roots equal to g x . Hence 
the axes corresponding to latent roots occurring only once among the latent 
roots of ty are never indeterminate. 

Identical Equation. 

§3. If <£> is a matrix all of whose latent roots are distinct, any quantity p 
in the extension of the ground may be represented linearly in terms of its axes. 
Thus for any quantity p we may put 

p = %pi + »2p? + • • • • + 2 M p«, • 

Operating upon p by <?> — g M we get 

(<?> — gj) p = % (ffi — sO pi + *% {g% — g.) p 2 + — + ^{g, — gJ)o m . 

By this operation the last component of p, that along p„, is annulled. Operating 
on (4> — g M ) p by <£> — g ia _ 1 , another component of p is annulled. Finally 

(* — ■&)(* — 9%) (<?> — fl r — 1)(* — 0.) P = ° • 

This being true for any quantity p in the ground, may be written 
(f—9i)(<l>—9%) (<?>— 9»-i)(<i> — &>)=<?> M — m«.-i4>" -1 + =F»ni^±«i= Q. 

This is Cayley's "identical equation." 

The proof of the identical equation may be extended to any case as follows : 
The latent roots being supposed to be all distinct, let \==g 2 — g lt h 3 = g 3 — g lt 
etc., h n ^z g n ^-g x , and let 4> — g x be denoted by ^ ; finally, let # denote the 
product (<P^-g n +i)(<p — g n +z) • - ' ' ("P — g*)> — the identical equation then be- 
comes 

[^ - Xh,^- 1 + etc.] x = 4"jt - Xh^-'z + etc - = Q • 

This equation is true however small the A's may be, provided the latent roots 
remain distinct. Consequently, as the A's diminish without limit, the first term 
of the second member i^Z mus t ultimately diminish without limit, and in the 
limit the identical equation becomes 

+ n Z = ti>-9iY&-9n+i) ■■■■ (*-&.) = <>. 
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Thus it is proved that the identical equation still subsists when any group of n 
latent roots becomes equal, the other latent roots being distinct. In the same 
way it may be shown that if any other set of the latent roots become equal, and 
finally whatever the relation between the latent roots, the identical equation 
subsists. 

Oayley has remarked that the identical equation may be represented as 
follows : 



«11 — <?> «12 « • • • «lo> 

a 21 ®22 <?> • • • • #2o> 



= 0. 



Or, as stated by Cay ley, the determinant of the matrix less the matrix itself 
considered as involving the matrix unity is zero. This relation Cayley denotes 
symbolically by Det (l<p— $1) = 0, where $1 signifies that <?> is to be treated 
as a scalar. I propose to denote <p considered as involving the matrix unity 
by $, when, with the notation previously employed for the determinant of a 
matrix, the identical equation may be represented by | <|> — |> | = . 

Corresponding to any latent root g x of 4> can, by §1, always be found an axis 
of <£; and upon this the effect of any operator <£> — g 2 for a latent root g % distinct 
from g x is by §2 only to multiply it by a scalar constant. Consequently no 
product of factors <?> — g, not containing q> — g lt can annul the axis p lt and hence 
no such product can vanish. Similarly for every other latent root. Whence the 
identical equation must contain a factor <£> — g for each distinct latent root. On 
the other hand, it was shown in §2 that when the latent roots are not all distinct, 
there may be more than one axis corresponding to a latent root g x which 
occurs more than once. In this case <£> — g x will annul an extension of order 
greater than unity, and the order of the identical equation will be lowered. 

It is obvious that there can be but one identical equation, and that if <p is 
subject to any other equation involving only the matrix ty , this must contain the 
equation 1 4> — ^ | = as a factor. 

Converse of <£>. 

§4. The latent function | ^ — g \ of ^ is obviously identical with the latent 
function of 4>. Consequently the latent roots of <£> and ^ are identical, and the 
identical equation in <jj> is 
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The axes of <p will be denoted by p lf p 2 , . . . . p w . They do not necessarily 
coincide with those of q>. Indeed, I will show later that the axis p ± of cf corre- 
sponding respectively to the latent root g x is normal to all the axes of <£>, 
pai p3» • • • p« corresponding to the remaining latent roots. 

On the Form Representing a Matrix. 

§5. If the latent roots of <£> are assigned, q> is subject to a conditions. If in 
addition a linearly independent axes are assigned to which the assigned latent 
roots are to correspond, <p is by definition completely determined. When the 
latent roots are given as not all distinct, and corresponding to them a axes, of 
which those corresponding to a set of equal latent roots are not in every case 
linearly independent, <£> is not completely determined ; an infinity of matrices 
can be found which will have these latent root and corresponding axes, and the 
same identical equation. 

To find the form or array representing the matrix <p whose latent roots 
ffn 9%, et c-) and corresponding axes p lf p 2 , etc., are assigned, we may proceed as 
follows. Let 

Pi = ^n«i + *i2«2 + •••• + x u a m , 



p u = aw*i + a^a.; + ....+ x ma a a . 
Denoting by X the matrix formed from the array of the x's, evidently 

p x = Xa x , Qz = Xa % , p 3 = Xa 3 , etc., 
and by definition 

<P OZixi) = gxXa x , $ (Xa z ) — g 2 Xa 3 , $ (Xa 3 ) = g^Xa^ , etc. , 
i. e. ($lXla u a 2 , a w ) = (X\ g x .... Ja lt a 2 , a B ) 

o g z .... o 



o o 



g« 



As this equation is true for each of the a's, it is true for any expression linear in 
them, and so for any quantity in the ground. Hence 

<t> = (x\ 9l o ....o ix- 1 ) 

g 2 







g« 
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This method is substantially identical with that given in Clifford's "Fragment 
on Matrices" by Oayley as a correction of Clifford's result. 

If a set of o linearly independent quantities other than the a's are chosen 
to represent the ground, the linear unit operator or matrix <p will, in general, be 
represented by another form or array. Two exceptions are, however, to be 
noted, the matrix unity and the matrix zero always have the same representa- 
tion. By §2, if the latent roots of ty are all distinct, the axes of q> are capable 
of representing the ground. If they are chosen for this purpose, the form 
representing <£ becomes 

(9i .... ) 
g % 







9» 



Evidently this form of a matrix gives an immediate proof of the identical equa- 
tion. It has just been shown how from this form to pass to that form of q> when 
the a's are chosen to represent the ground, the axes being given in terms of the 
a's. If the a's are given in terms of the axes, and a 1 = Xp 1 , a % — Xfy, etc., the 
only change in the formula for <£ is that JTis to be replaced by X~\ 

The more general problem, given the form representing a matrix q> in terms 
of certain a linearly independent quantities (i. e. the form of <£>, when these 
quantities are chosen to represent the ground), to find the form of q> in terms of 
another set of a linearly independent quantities is solved in a similar way. 



Denote by <|> a the form of 4> in terms of the set (a 1( a 2 
form of <?> in terms of the set (/?!, /3 2 , . . . . /? w ), and let 

<pa — ( «ii <h% etc - ) % = ( &n 
etc. 



a M ) , and by fy, the 



&12 
^21 ^22 

etc. 



etc. ) 
etc 



If 



& = ajuttx + x n a z +.... + a? la) a„, 
& = <%«! + *22«2 + + aw* M , 



Ao = «W<*i + aw*s + •••• + « WM a w , 
and X represents the matrix formed from the array of x's ; then 
/3i = Xa 1( /2 2 = JTa 2 , (3 s = Xa 3 , etc. 
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Operating by the matrix ^> upon the /?'s and their equivalents gives 
faX^^fa^x, faXa % — fa^, faXa 3 = fa(3 s , etc. 

For, by definition, the effect of fa upon any expression in the a's is identical with 
the effect of ^ upon the equivalent expression in the /?'s, since fa and fa both 
represent the same matrix. It is not proper yet to replace the (3's by their 
expressions in terms of the a's, for fa as an operator upon the a's is not equiva- 
lent to the effect of the matrix <£>. As the result of operating upon the /3's by 
fa we have 

faXa,! = Jn/?! + 6 2] /? 2 + ....+ 6 wl /3 M , 

faXa % = 5 12 /?i + b n (3 % + .... + & w2 /? w , 



faXa a — b la ^ + & 2w /3 2 + ....+ b^Pn. 
Replacing now the /3's by their equivalents, 

<?> a Xa x = 6 n Xax + & 2 i-Xa 2 + ....+ & M i-2a«> = Xfaax, 
faXa% = b n Xai + & 22 Xa 2 + ....+ b a% Xa a = Xfaa%, 
etc. etc. 

i. e. ^aJTCai, a,, . . . . a„) = Xfa(a 1; a 2 , . . . . a„), 

where ^ is now regarded as an operator upon the a's, 

.". fa = J^jJf -1 . 

§6. From the definition of a matrix as an operator upon a certain ground it 
follows that the identical equation does not vary with the form. It may, how- 
ever, formally be proved as follows : Since fa = XfaX -1 , evidently $• = XqfyX~\ 
etc>, let the identical equation in fa be 

K — »»— lflT 1 * ^=m 1 fa±m =0, 

.-. X(<^ — m^!^" 1 + =f m^ ± m) Xr 1 = 0. 

But since the /3's constitute a set of w linearly independent quantities, \X\ =f=0. 

Consequently, 

$? — m a _ 1 fa- 1 + =fm 1 fa±m = 0. 

Whence the vacuity of a matrix does not vary with its form. This proof of the 

in variance of the vacuity does not, however, take into consideration the fact that 

the identical equation may degrade. The invariance may, however, be shown as 

follows : If g is a scalar XgX~ 1 = g ; hence 

fa - g = XfaX- 1 - XgX- 1 = X(fa - g) X~\ 
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From the definition of the product of two matrices it follows that the determi- 
nant of their product is equal to the product of their determinants. 

.-. \<p* — g\ = \X\ \^ — g\ [X- 1 ] = \% — g\; 

for | X -1 1 = | X\ -1 . Hence not only are the distinct latent roots of ty unchanged 
with its form, but the number of times each latent root is repeated is also 
unchanged. 

Definition of Nullity and the Law of Nullity. 

§7. A null matrix is one whose determinant vanishes, or of which all the 
minors of a certain order vanish. A non-null matrix is said to have a nullity 
zero, and one, every constituent of which is zero, is said to be absolutely null, or 
to have the nullity a. It has been shown that the absolutely null matrix is the 
scalar quantity zero. Between these limits the number expressing the measure 
of nullity may have any integer value. If all the (x — l) th minors of the deter- 
minant of a matrix vanish, but not all the x th minors, the matrix has a nullity x . 
Nullity of order one or simple nullity is evidently the same as simple vacuity. 
The vacuity of a matrix obviously cannot exceed its nullity, but it may have 
simple nullity and vacuity of any order from unity to a. 

The nullity of 4> is not affected by multiplying it by a non-null matrix. 
Thus if the nullity of <£> is «i, the nullity of & zero, the nullity of <£>» = 4* is m. 
For the (m — l) th minors of 4> consist of all possible products of a rectangular deter- 
minant formed from a — m + 1 rows of $ into the rectangular determinant 
formed from the corresponding a — m + 1 columns of a ; and each of these 
products is resolvable into the sum of a — ra + 1 products of an (m — l) th minor 
of 4> into an (m — l) th minor of m. But the (m — l) th minors of q> are all zero. 
Whence the nullity of 4- is not less than m . 

Since, however, | w\ =j= 0, q> =^m~ 1 , and in the same way it may be shown 
that the nullity of <£> is not less than the nullity of 4 1 ; hence the nullity of 4 1 is m • 
By the same method it may be shown that if m$ =,4', the nullity of ^ is still 
equal to that of $.* It follows immediately from this that the nullity of a 
matrix is unchanged when the form representing it is changed. 

* This proof follows the method employed by Sylvester to prove the theorem regarding the lower 
limit of the nullity of the product of two matrices. To Sylvester, as has been stated (I, §10), the term 
nullity and the discrimination between degrees of nullity are due. 
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§8. If the greatest extension annulled by any matrix is ni, its nullity is m, and 
conversely. This I call the first branch of the law of nullity. Let a x , a 2 , . . . . a m be 
any m linearly independent quantities in the greatest extension A annulled by <£>; 
and let /? X) /? 2 , ..../?„ be any n= a — m quantities in the extension B comple- 
mentary to A with respect to the ground ; the a's and (3's may then be chosen to 
represent the ground, and any quantity o in the ground may be expressed by 

+ aw*. + yiPi + yA + ....+ y n {3 n , 

+ t/n$Pn- 

+ a ml a m + bnp! + b 2J fa + + b nl (3 n , 

+ ct^a m + b n $ x + hA + + &*«/?», 



If 



■• 4>p = yiPi + y&fo + 

,<p{3 1 = a u a x + a 21 a 2 + 
<p(3 2 = a 12 a x + « 22 a 2 + 



<p(3 n = a^tti + a 2w a 2 + . . . 
then 4> is represented by the form 



+ a mn a m + »i M pi+ t»2«P2 +.... + b n J3 n , 





a u 

«21 


a lz . . . . a ln 
o 22 . . . . a 2 „ 


a «ll 


^m2 * • ' * ®mn 




^21 


6 12 .... 6 lB 
6 33 .... 6 2n 




Ki 


"»2 .... O nm 



where the blank square and rectangle consist of wi columns of zeros. Hence <p 
has a nullity at least ni . If the nullity exceeds m , the determinant A formed 
by the n columns of the 6's must vanish, and likewise every determinant which 
can be formed from n == a — m rows of the last n columns of <?>. Denoting the 
first minors of A, corresponding respectively to b u , b 12 , etc., by B n , B 12 , etc., we 
shall then have 

<?> (B li S 1 + B, A +....+ B u (3 n ) = Ai^i + BrfP* +••••+ B ln ^ n = . 

But by supposition the extension A contains all the axes of ty corresponding to 

47 
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the latent root zero.* Hence the nullity of 4> cannot exeeed m. The converse 
is immediately evident. 

§9. If the matrix <f> transforms every quantity in the ground into one in an 
extension of order a — in, the nullity of q> is m, and conversely. For if 

p = x x a x + x 2 a 2 + '....+ xjt M 

is any quantity in the ground, and 

<?>P = (^n^i + «u«8 + + <hJX«) <*i + (an x i + a n x 2 + ....+ a u xj) a % 

+ .... + (a^ + a^x z + ....+ a lm a; w ) a M 

is contained in an extension of order a — m, then there are m linear relations 
between the coefficients of the a's in the last equation true for all values of the 
as's. These give m sets of a equations between the constituents of the matrix ; 

namel y - ?i«n + w + .... + ki = o , 

'1^12 T" '2^22 T" • • • • "r 'iAo2 ~~ "» 
n a u "T 4 a 2o T" • • • • T IjX'wo — - > 

and m — 1 other sets of a equations, in which the Vs are replaced in each set by 
other constants. Prom the w equations obtained by taking the first equation in 
each set, m — 1 of the a's may be eliminated, leaving ,one equation between any 
a — m + 1 of the constituents of the first column; and, similarly, an equation 
linear in the corresponding constituents of the second column may be obtained 
from the second equations of the m sets, etc. Thus result a equations linear in 
the constituents of each column taken from any a — m + 1 rows of | $ | , the 
coefficients of the a's being the same for all constituents in the same row. The 
resultants of all sets of w — m + 1 of these equations, together with the like 
resultants formed from all such groups of a equations, constitute the (wi — l) th 
minors of \q>\ ; and these resultants must be zero. Whence the proposition. 
The converse is evident. 

If, as in the preceding section, A denote the total extension annulled by <p , 
and B the complementary extension, it will not in general happen that B is the 

* If not only A but also all the JB's vanish, the analysis needs correction. But, since the constituents 
in the last n columns of <p do not all vanish, the minors of order p {for some value p from p = 1 to 
jj = w — 1) do not all vanish ; and then in a similar way, if the minors of order greater thanp all vanish, 
it may be shown that there would result a linear relation between ^j , p/3 2 , etc., i. e. <!> would annul a 
certain expression linear in the /3's. 
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extension into which $ transforms any quantity of the ground. As will appear 
later, this is only true when the vacuity of <p does not exceed its nullity. Let B ± 
denote that part of B transformed by $ into A. There must exist such an 
extension unless no extension but A is annulled by any power of <p. For if <p K+y 
(where x =1= 0) is the lowest power of <?> that annuls an extension G other than A , 
then <p*G is included in A, since by supposition A is the total extension annulled 
by <£> ; moreover, ^» K ~ 1 G is not included in A, since it is not annulled by <£>; hence 
^> K ~ 1 G is a part of the ground complementary to A that is transformed by 4> into 
A. Let B 2 denote that part of the extension B — B 1 which is transformed by ^> 
into B 1 . By an argument similar to that above, it may be shown that if the 
order of B % is zero, then no power of 4> annuls an extension other than the 
aggregate of A and B{, i. e. A + B x ; etc. Finally, let B p denote that part of the 
extension B — (B x + B % + ....+ B p _ x ) which is transformed by <£> into 5 J> _ 1 ; 
and suppose no portion of the remaining extension 

B p + , = B - ( B, + B z + . . . . + B, _ x + B p ) 

is transformed by q> into B p . Then all the extensions A, B lt I? 2 , .... B p _ lt B p 
are annulled by some power of <J>, but no part of the extension B p + 1 comple- 
mentary to their aggregate is annulled by any power of $. The following mul- 
tiplication table shows the effect of <£> upon these mutually exclusive extensions: 



A 



B, 



B, 



etc. B, 



-i 



B n B 



'p+i 



¥ 



etc. 



f 



,*>-! 



4>p 






A' 


B{ 


.... 


*',-> 


SJ-i 


B P+ i 








A" 


.... 


B'J- S 


BU 


*i'+i 











« • * • 


Dill 
■°p-i 


-Dill 


B'JU 


etc. 

























JlP-1) 


Bff-v 


B&? 
















A™ 


mi 
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In this table the accented letters denote extensions included each in the exten- 
sions denoted by the same letter with fewer accents or unaecented. 

The total extension A of order m annulled by q> I shall term the null exten- 
sion or null region of the matrix, and the aggregate of the extensions exclusive 
of A annulled by some power of q> the vacuous extension or vacuous region of the 
matrix. The extension or region into which <£> project any quantity of the 
ground is 

Q(A + B 1 + B t +.... + B p _, + B p + B p+l ) 

= A> + B> + B' + .... + B p _ 1 + B p+1 ; 

and this I shall term the residual region of $ . If { A 1 } denote the order of the 
extension A', etc., evidently 

{*}<&}, \B{}<\B % \ etc. {B'^S^l \B> + 1 \<{B p + 1 \* 
p—i p+i 

.-. {A'\ + J2 {B> K \ + \B p + 1 \<J2W=»- 



m. 



But by the proposition just proved the order of the residual extension is a — m, 
.: \A'\ = {B 1 \, \B[\ = \B 2 }, etc {B P _ 1 } = {B P \, {B p + 1 \ = \B P+1 }, 

and so B p+1 = B p + x . 

Hence if G is any extension which has no part in common toith the null region 
of 4> , it is transformed by 4> into an extension of the same order. 

§10. The second branch of the law of nullity is easily derived from the two 
preceding sections. Denoting as before by A the null r-egion of ^ of nullity m, 
and by B the extension complementary to A with respect to the ground, let 
denote the null-region of ^ of nullity n, and let |p denote that part of B which 
4> transforms into O. Obviously, the null-region of 4«p is the aggregate of A, 
the region which (p annuls, and of |p , that part of the ground complementary to A 
which <|> transforms into G. By what was proved in §9 it follows that the order 
of ^)^ is the same as that of || . Whence to determine the order of g it suffices 
to determine the order of 4>|p , the extension common to G and the residual 
region of <|>. The order of the residual region of ty is a — m and the order 
of G is n; hence the extension 4>|| common to G and the residual region of q> 

* It is obvious from the definition of a linear unit operator that it cannot increase the order of an 
extension ; so if <t>A ~ B , hence \B\S { A \ . 
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is at most of order n. If (a — m) -\-n ~^>a, i. e. m<C,n, the extension G and the 
residual region of q> have common an extension of order at least ((o — m) -\- n) 
— o = n — m, but if (g> — m) + n<o>, i. e. m>n, these two extensions do 
not necessarily have any part common. Hence the order of q>^, and conse- 
quently of |p, is at most n; and if m<in, the order of §} is not less than 
n — in, but if m>n, the order of |p may be zero. The extension A is of 
order m and has no part in common with ||. Hence the null-region of 4rf> 
cannot be of order greater than m -\- n; and if «i<C n, it cannot be less than 
(n — rn)-\-rn=n, but if m>n, it cannot be less than in; and, consequently, 
the nullity of the product of two matrices is not greater than the sum of their 
nullities nor less than the greater nullity of the two matrices* This theorem is 
due to Prof. Sylvester, who terms it the law of nullity. Owing to the impor- 
tance of the relation of the null-extension to the nullity of a matrix, I term the 
whole relation of the null-extension and nullity of one or more matrices the law of 
nullity, and this theorem the second branch of the law. 

Suppose that G, the null-region of 4s has no part in common with the 
vacuous region of <p; let E denote the extension common to G and A, the null 
region of <?>; then the residue of G, namely G — E, is wholly contained in that 
part of the ground which is complementary to the aggregate of the null-region 
of <£> and the vacuous region of <?>, and which was denoted in §9 by B p + 1 . Let 
Hi denote that part of the ground complementary to A which is transformed by 
^> into E (Hi is evidently contained in B x ) ; let the order of E be p and the order 
of Ipx be q. Since <pB p + 1 = B p + 1 (§9), the null-region of •$>$> is obviously the 
aggregate A, of ^, and of G — E, which are all mutually exclusive ; and, con- 
sequently the order of the null-region of 4"p is m, + g -f- (n — p). Hut p<m, 
p<n; moreover, the order of the extension tp^ is equal to that of ||i, but Q^ 
is included in E, hence q^p- If p = , the order of null-region of 4«p is m + n . 
In this case the null-region of 4* is included wholly in that part of the ground 
complementary to the aggregate of the null-region of q> and the vacuous region 
of i^. The aggregate of these two extensions I shall term, for reasons which 

* The outline of this proof was communicated by me to the Johns Hopkins University Mathematical 
Society in a paper read before the Society, Nov. 1888, when I was informed by Dr. Franklin that this method 
of proving Sylvester's law had been employed by Mr. Bucheim in a note in the Phil. Mag. (XVIII, 459). I 
had previously had no knowledge of any anticipation of this method of dealing with the theory of matrices 
except what is contained in Clifford's " Fragment on Matrices. " Subsequently I examined Mr. Bucheim 's 
proof, and I find that though it is sufficient to give the limits between which the nullity of the product of 
two matrices lies, nevertheless it is defective in assuming what is not always true, viz. that the resultant 
region of a matrix is coincident with the extension which in §9 I have denoted by i? p+ i . Mr. Bucheim 
also omits to prove that the order of 18 is the same as that of <*33 . 
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will appear presently, the latent region of q> corresponding to its latent root zero. 
Hence if the null-region of 4 1 has no part in common with the latent region of ^> cor- 
responding to the latent root zero, the nullity of 4^> is the sum of the nullities of 
<|> and 4'. A fortiori, if the latent regions of ^ and ^ corresponding to the latent 
root zero are mutually exclusive, the nullity of 4«P (and also of q4>) is the sum of 
the nullities of <£> and 4 • 

If the order of the vacuous region of <J> is zero, so that the null-region of q> 
is coextensive with the latent region of <J> corresponding to the latent root zero, 
the null-region of 4^ is the aggregate of A and 0; hence if these extensions do 
not intersect, the nullity of i'Q is m + n. Conversely, */ the order of the 
vacuous region of 4> is zero, and the nullity of <fyty is in + n, then ilie null-regions of 
q> and 4 1 have no part in common. 

Nullity of the Factors of the Identical Equation. 

§11. If the latent roots of q> are all distinct, the nullity of the product 
4>j = (<£ — g^)(<p — (ft) . • . . ($ — g m ) of m factors of the identical equation is m.* 
This theorem is Prof. Sylvester's, and is termed by him the corollary of the law 
of nullity. His demonstration is as follows : Since each factor, being simply 
vacuous, has a nullity of order unity, the nullity of the product $> x cannot exceed 
m . Similarly the nullity of the product <J) 2 of the remaining factors of the iden- 
tical equation cannot exceed a — m . But the nullity of the product of c&j and <J> S 
is a, and consequently the sum of their nullities must be as great as a. Hence 
the nullity of 4> 2 cannot fall short of m , and the nullity of <E> 2 cannot fall short of 
a — ni. 

When the latent roots of q> are not all distinct, the law is not so simple. 
Suppose the distinct latent roots of q> be * in number, namely, g x , g 2 , . . . . g t , 
occurring severally m 1 , m % , . . . . m t times, and that the identical equation is 

(<?> - 9i) mi ~ Kl (<£ - 9%)"'" • • • • (<?> — 9^"* = 0. 

For the investigation of this case the following lemma is required, namely, that 
the vacuity of any positive integer power of <£> — g, for any latent root g, is equal 
to the vacuity of $ — g, and hence to the number of times the latent root g 
occurs. This lemma is an immediate consequence of the theorem that the 
determinant of the product of two matrices is the product of their determinants; 

*From the law of latency, §15, it follows that the vacuity of * t is also m. 
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for then, if m is a positive integer, if (i is a primitive 171 th root of unity, and if 
^ = <?> — 9, 

\V~g*\==\+-g\.\+-i t g\.\+- l fg\....\4,- ll ~-*g\. 

But if g r is the lowest power of g that appears in the latent function | ^ — g\ of 
4s then (ng) r is the lowegt power of (ig that appears in \4> — (ig\, etc.; whence 
the lowest power of g m that appears in 1 4> m — g m \ is the r th power of g m . Conse- 
quently 4 ,m = (^> — gi) m has the same vacuity as $ = <p — g x .* Resuming the 
investigation of the nullity of any product of matrices 4> — g for different latent 
root g, as a consequence of this lemma, the several factors (<|> — g^) mi ~ Kl , 
(q> — g%) mi ~ H , etc., of the identical equation have the, vacuities m 1 , m z , etc.; and 
hence their several nullities cannot exceed vn x , #? 2 , etc., respectively. But the 
nullity of their product, which cannot exceed the sum of their nullities, is a. 
Hence the nullities of the several factors are n? 1} m z , etc., respectively; and 
hence corresponding to each latent root g of (p is an extension of order equal to 
the number of times that latent root is repeated which is annulled by that power 
of $ — g appearing in the identical equation. f These extensions I term the 
latent extensions or latent regions of the latent roots. In a similar way it may 
be shown that the nullity of the product of any two factors of the identical 
equation (<p — gi) mi ~" 1 and (<£> — g^ m *~""' is m 1 -\- m 2 ; but since the vacuity of 
either matrix is equal to its nullity, the order of the vacuous region of either 
matrix is zero; and hence by §10 their null-regions have no part in common. 
Similarly with respect to the null-regions of any two other factors of the 
identical equation. Hence the latent regions are mutually exclusive.^ Since 
each positive integer power of 4> — g x is vacuous, there is, evidently, an exten- 
sion of order at least unity annulled by ty — g x ; an extension annulled by 
(<p : — g x y, etc. ; and each of these extensions is included in the one corresponding 
to the next higher power of 4> — g 1 in the series, since if an extension is annulled 
by any power of <£> — g 1 it is annulled by the next higher power. But as the 
nullity of no power of <p — g x is greater than m x , which is the order of the latent 
region corresponding to g lt all these extensions are included in this latent 
region. Similarly with respect to the other latent roots. The latent region of 

* It should be remembered that the latent roots of a matrix are the roots of its latent function. 

t In this method of finding the nullity of the several factors (f — g s ) mi ~ " 1 of the identical equation, 
I have followed another one of Sylvester's methods of demonstrating the corollary of the law of nullity 
when the latent roots are all distinct, the only case Sylvester considers. 

J That the latent extensions' are mutually exclusive may be proved very simply by the method 
employed in the first part of ?12. 
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q> corresponding to the latent roots g x ,g z , etc., are, however, obviously the respective 
latent regions corresponding to the latent root zero of the matrices 4> — g lf $ — g % , 
etc.; and consequently these extensions are also respectively the latent regions 
corresponding to the latent root zero of any positive integer powers of <£ — g x , 
q> — g 2 , etc. Whence the nullity of the product whose factors are powers of 
q> — g x , <£> — g 2 , etc., is by §10 the sum of the nullitiesof the several factors, as 
the latent regions of these factors corresponding to the latent root zero are 
mutually exclusive. 

The next problem is to find the nullity of successive powers of q> less than 
any of its latent roots ; and I shall show that the nullity of the (<p — g x ) 2 is greater 
than the nullity of q> — g 1 by an amount at least unity (unless the nullity of 
Q — g l is m ± ) , and that the nullity of successive powers of ^> — g x goes on increas- 
ing by an amount not greater than the increment of nullity in the preceding power, 
until some power of <£> — g x is reached whose nullity is equal to its vacuity ; that 
power of 4> — g x whose vacuity is equal to its nullity is evidently the factor in 
the identical equation corresponding to the latent root g x . This theorem has 
already been proved in §9, where it was shown that the region annulled by 
(4> — g x f consisted of A the null region of $ — g lt and of B x the extension trans- 
formed by 4> — g x into A, of order equal to or less than that of A, etc. The 
theorem may also be proved as follows : Denote <p — g 1 by $, let the null region 
of 4* of order a be denoted by A ; and, of the complementary extension, let B x 
of order b be that part transformed by $ into A. Such an extension must 
exist, otherwise 4* would annul only the extension A; hence A would also be 
the null region of 4 /3 , since 4? can annul only the null region of 4? together with 
that extension projected into it by 4', etc. Finally, no power of 4> would annul 
an extension other than A ; but then the factor -v|/ would occur in the identical 
equation only to the first power; consequently the nullity of ^ would be as great 
as its vacuity. The order of B x cannot be greater than the order of A , for then 
the ^ of any b linearly independent quantities (3 1} (3 t , . . . . (3 b , in the region B x 
would be expressible linearly in terms of a <C b linearly independent quantities 
in the region A ; hence for some value of the t's other than all zero, the expres- 
sion ^ + ^ +.... + t^ b = ^ («,& + t£, + . . . . + b ) 
would vanish, which is contrary to supposition, since no part of B x is in the null 
region of i^- In the same way it may be shown that if the nullity of 4 ,m+1 is not 
greater than the nullity of ^ m , no higher power of 4> has a nullity greater than 
the nullity of 4 ,m > which must therefore be equal to its vacuity; and also that the 
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extension annulled by 4' n + i additional to that annulled by ^ n+1 cannot be of 
order greater than the order of the extension annulled by $ n+1 additional to the 
null region of 4 ,n - 

§12. No rational integral function of <p operating upon any quantity in the 
latent region corresponding to one latent root can transform it wholly or in 
part into the latent region of any other latent root. By Sylvester's formula, §16, 
it may be shown that any function of $ may be reduced to a rational integral 
function ; whence the above proposition may be stated for any function of <£> . It 
is obviously only necessary to prove that no function of <£> can transform any 
quantity in one latent region wholly into, the latent region of another latent root. 
Suppose F<p is a rational integral function of 4> of order n ; let £ x be a quantity 
in the region of g x . If F$.i$ 1 = £ 2 , where % s is -a quantity in the latent region 
of g % , then for some integer rn>tn z — x 8 , 

= (!p-9t) m &=0. 

Suppose F<p contains <£> — g 1 as a factor to the <p tb power, then unless 
F<p%x = , a case which need not be considered since £ 2 is to be regarded as non- 
evanescent, there exists a linear relation 
(4> -&)»+"&+ 4 ($~- fl r 1 )"+"- 1 &+ . . . . +L(<p-9iY+% + M(<l>-g i n 1 = 0, 

where M=j=0, between ($ — <fr) p £i and other quantities that are annulled by 
successive powers of 4> — g x , which is impossible. 

Since the latent regions are mutually exclusive and together make up the 
extension of the ground, a number of linearly independent quantities may be 
taken from each latent region and together may be employed to represent the 
ground. 

By what has just been proved, each latent region in respect to <£> constitutes 
a subordinate ground ; for the effect of the matrix $ upon any one of these 
extensions is to convert a set of linearly independent quantities, equal in number 
to its order, into other quantities in the same extension. It would thus seem 
that 4> might be regarded as an aggregate of subordinate matrices corresponding 
to and equal in number to the distinct latent roots, each of which would have 
upon the latent region corresponding to it the effect of ty , and would have a null 
effect upon any other extension. This suggestion may be verified as follows : Let 
the latent roots of $ be i in number, namely, g lf g % . . . . g t , occurring severally 
tn, n, . . . p times ; let the latent regions corresponding to the g's be A, B, ... L 
48 
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respectively ; and, finally, let any m. linearly independent quantities a 1( a 2 , . . . a m 

be selected from A, any n linearly independent quantities (3 lt (3 Z /?„ be 

selected from B, etc., and let any p linearly independent quantities \, h z , . . . \ 

be selected from L; these quantities together embrace the extension of the 

ground. If 

<£> ai = a llttl -f- a n a % + . 

4><*2 = «12«1 + «22«2 + • 

<P& = haPi + &J2& + • 

GlC.j GcC.j 

^i = ^11^1 + h&% + • 

<^g = Zjg^j .+ ^2^2 "h • 

then the form of $ becomes 



T" ^m2 a mi 6tC, 
"T OnlPnt 

+ 2W*„, etc., 
i vpiA/p i etc., 



Ct-11 d\% .... Cli m 
^21 ^22 • • • • ^2m 








^«ll ^m2 « • « « ™mm 




b u b n . . . . 6 ln 

2 l 0-32 • •'• • #2» 






Oril ®n1 • • • ' ®nn 
















hi In ... . hp 

hi '22 • • • • 4p 


*pl (pj . . . . 6pp 
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In this form all the constituents except those in the squares along the diagonal 
are zero. It is obvious from this form that the array of the a's forms a matrix 
by itself, and that its effect upon any region other than that of the latent root 
g x is to annul it, etc. Whence if fa denote the matrix formed from the a's, and 
fa that formed from the 6's, etc., then fa, fa, etc., are nil-factorial together, and 



and, in general, 



<?> = $1 + #2 + • • ■ • +fa, 

^ = # + #+....+#; 
Up = Ffa + Ffa + + Ffa. 



The matrices fa, fa, etc., regarded as pertaining to the ground of <£>, are, of 
course, vacuous ; thus, trom this point of view, fa has the representation 
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a n • • 


• • ^lm 


o . . 


..0) 


a 21 


a 22 . . 


• • ^2m 


o . . 


. . 


U>m\ 


a m % . . 


• * ™"mm 


o . . 


. . 





. . 


. . 


o . . 


. . 





. . 


. . 


o . . 


. . 



and has a nullity at least co — m; if ^> is non-vacuous, the nullity of fa is just 
a — rn, etc. But regarded as subordinate matrices of different systems, per- 
taining to the subordinate grounds A, B; . . . . L, then, unless <£> is vacuous, 
fa, fa, .... fa must be considered as non-vacuous matrices; and fa will have 
the representation 



On 

«21 



a n 



a 






)> 



^ml drnZ 



a. 



'mm A 



this matrix being supposed to operate only upon that part of the ground com- 
prised in the subordinate ground J., namely, the extension of the set (a x ,a % ,... a m ) , 
and must be regarded as not susceptible of operating upon any other ground. 
The subscript A is employed to denote that the ultimate operands of this matrix 
are only expressions linear in the a's, etc. From this point of view it is proper 
to consider $!, fa, etc., as having reciprocals, providing 4> is not vacuous; and 

then ^- 1 = 4»r 1 + ^ 1 + .... + <?>r 1 . 
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If <j!> is vacuous it must have zero as a latent root ; but only one of the 
subordinate matrices corresponds to this root ; consequently one and only one of 
the subordinate matrices will be vacuous, and evidently the vacuous subordinate 
matrix will have exactly the same vacuity and nullity as q> . The vacuous subor- 
dinate matrix is evidently a nilpotent quantity, since if <p is raised to a suffi- 
ciently higher power it annuls the region corresponding to the latent root zero. 
Thus if <7x= 0, then A is annulled by ty™ 1 *" 1 ; hence 

f--"' = o + fy™ 1 -" 1 + ....+ $r i-Ki - 

The conception of a matrix as a sum of a set of subordinate matrices is more 
readily grasped by the consideration of Peirce's linear-form representation 
of a matrix. Thus 

4> = 2£ a a rs (a r :a,) + 2 r 2 8 b r8 (/V-&) + + 2£ 8 l r8 (K '%>*)• 

But evidently the set of vids (a r : a 8 ) form a quadrate system by themselves, simi- 
larly with respect to the set of vids (/3 r :/? s ), etc., and any linear expressions in 
the vids (a r :a s ) will have as ultimate operands only expressions linear in the 
a's ; while any expression linear in two different sets of vids, if they are to be 
regarded as susceptible of operating upon each other, are mutually nil-facient and 
nil-faciend. 

Considering the vids (a^:a g ) as a quadrate system by themselves, scalar unity 
will be expressed by 2fr (a^:a r ) , and may be denoted by l x . The quadrate system 
formed from the vids (/3 r : /?„) will have as its scalar unity 2"r (fi r : /?,.) , which may 
be denoted by l s , etc. We evidently have 

Ms = Mi = lil 3 = Mi = etc. = ; 

and, denoting by 1 the unity of the complete system, 

l = l 1 +l,+ .... + l i . 
Thus — g = (ft — glj + (fr — flrl,) + ....+($«- gl t ). 

Since any matrix can thus be resolved into an aggregate of as many other 
matrices (mutual nil-factorial), as it has latent roots, each subordinate matrix 
corresponding to a latent root and being of order equal to the number of 
times that latent root occurs, hence in general it suffices to prove a theorem 
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relating to a single matrix for one all of whose latent roots are equal, when, if 
true, it may by means of this proposition be inferred for any case.* 

§13. Having obtained the law governing the nullity Of the factors of the 
identical equation, it is now possible to solve the problem touched upon in §3 
with regard to a matrix subject to a condition involving only itself. As was 
stated above, it follows from Sylvester's formula that any condition equation to 
which a matrix is subject may be reduced to a rational integral equation. Let 

Fp=($-gtf($-g t y. . . . ($-g m y. . . • ($-g n y=0 

be the rational integral equation expressing the condition to which <£> is subject ; 
then <p will satisfy this condition if any vn < w of the gr's are selected as its latent 

roots, as g lt g % g m , each occurring a', (3', .... a' times respectively; and if 

(<?> — giY' - " 1 nas tlie nullity a 1 , (# — g^Y'-"* the nullity /?', etc., and (<?> — g m Y'~" a 
has the nullity a', provided the sum of the accented Greek letters is a , and 
a' — %i < a , /?' — x 2 <(3, etc. For then 

(* - 9iY'- K > dp - mY~ K ° ••..(*- g-Y-~ = o , 
••• ^—9iY^-9»Y- • • • ($-g m Y- ■ • • ($-9nY=o. 

Conversely, if <p satisfies the last equation it must have a certain number of 
the gr's as latent roots, each occurring a sufficient number of times to make the 
total number of latent roots equal to a; and if g is any one of these latent roots, 
some power of <?> — g equal to or lower than that which appears in the above 
equation must have a nullity equal to the number of times g occurs. 

Law of Qongruity and Law of Latency. 

§14. If px is an axis of ty, it is an axis of any integer power <p m of <£; for if 
$p 1 = <7 1 p 1 , hence ty m oi = giQx . Similarly with respect to the other axes of $. 
Hence the axes of <p are all axes of q> m ; but since the axes of q> may be linearly 
related, they do not necessarily constitute all the axes of q> m . In like manner 

the axes of any root q> n (where n is an integer) are comprised among the axes 
of <£. We have for each latent root g r of <£>, v being a primitive n th root of unity, 
the equation 

* Thus if the identical equation has been proved for a matrix of any order whose latent roots are all 
equal, it may be shown by the principle of this section to hold in any case. For if 

= 01 + 02 + • • • • + 0il 

and = (0! — flnli)™ 1 = (02 - ff2l 2 ) m2 = etc., 

.-. (* -ft)- (0-gf 2 r* .... (#-0,)», = O. 
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(<P - 9?){V - *#?)(*" — v *9?) (<F - *> n ~ Vr) ?r = (* — fr) p, = . 

Evidently one at least of the matrices which in the first member operate upon 
p r must be vacuous, and hence one of the « th roots of each distinct latent roots 

of p is a latent root of p n . If then the latent roots g x , g % , etc., of p are all dis- 

tinct, so also are the latent roots of p n ; and then, by §2, the axes of p, p 1( p 2 , 

etc., are a linearly independent and determinate quantities, as are also the axes 

a 
of 4)""; and since the latter are comprised among the axes of p, it is evident 

that every axis of p is an axis of p n . Hence every axis of p is an axis 

m m m m m 

of <£"; and for the a axes we have the equations p»p 1 = gjp 1 , p™p % = gfp 2 , 
etc. It may now be established by the method of limits, in conformity with 
the definition of I, §8, that for any scalar m we have for the a axes the 
equations p m pi = gTpi, $ m p» = gsp»t etc. Consequently 2&4> m .pi = 2%f.p 1( 
2&4> m p 2 = 27tf<7™.p 2 , etc., provided the coefficients and exponents are scalars ; 
and so if F<p is any function of q> formed by the addition of scalar multiples 
of scalar powers of <£>, then for the a axes F$p 1 = Fg 1 .p 1 , F<pp 2 = Fg z .p z , etc. 
Thus what has been proved is, when the latent roots of q> are all distinct, every 
axis of <£> is an axis of Fq> ; and since the above equations, or their equivalent, 
(F<p — Fg^) pi = , (F<p — Fg 2 ) p a = , etc., hold only for the latent roots of Fp , 
hence tohen the latent roots of p are all distinct, the latent roots of Fp are the same 
function of the latent roots of p . 

The first of these theorems I term the law of congruity of the axes. The 
second is very important ; it is due to Prof. Sylvester, who terms it the law of 
latency. In the next section I shall extend the proof of the law of latency to the 
general case. The term law of congruity was employed by Prof. Sylvester 
interchangeably with the term law of latency. 

The converse of the law of congruity is also true, subject to a slight modifi- 
cation when two or more of the latent roots of Fp become equal. In this case, 
if the axes of Fp are first selected, they may not all prove to be axes of p ; but 
when this occurs, other quantities may be chosen which, together with the axes 
common to both p and Fp , shall constitute the a axes of Fp and also be axes 
of p. Thus, to take an example, let p lt p 2 , p 3 be the three axes of the ternary 
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matrix <£> whose corresponding latent roots are 1 , — 1,2 respectively ; then the 
axis of 4> 2 corresponding to the latent root unity is indefinite ; for 

¥ (xpi '+ ypt) = xpi + yo % . 

In general, this quantity will not be an axis of $. However, p 1( p 2 and p 3 are 
axes of both (jb and ^> 2 . 

If the latent roots of <p are all distinct, the ground may be represented in 
terms of the axes of $ ; in which case, if Ftp is any function of <p involving only 

the matrix 4> and unity, we have obviously 

F<p = (F gi ....0 ) 
Fg 9 ....0 



Fg„ 



From this form of F<p immediate proofs for the general case of the law of latency 
and of the law of congruity may be derived. 

§15. If the latent roots of <p are all distinct, it is not necessarily true that 
the latent roots of F<p are all distinct. When this is not the case, the order of 
the identical equation in F<p is lowered by unity for each latent root of F<p which 
becomes equal to another. Thus suppose Fg x = Fg z = .... = Fg n , then 
pi i pa> • • • • p n are a ll annulled by Fq> — Fg^ consequently the application of 
this operator to the most general expression p in the ground (which can now, by 
§2, be expressed in terms of the a axes of <£) annuls n components and leaves 
a — n components which will be annulled by the remaining factors of the iden- 
tical equation in F(p . I. e. the degree of the identical equation is lowered by 
n — 1. 

The law of latency may be proved in the general case as follows : Let the 
latent roots of ^> be all distinct, and let g % — g 1 =.h 2 , etc.; as h diminishes, 
Fg z = Fg x + \F'g x + etc. approaches Fg x . In the limit <£> has two equal roots 
g lt and Fq> has two latent roots equal to Fg v Obviously, when <£> has two or 
more equal latent roots, the latent roots of Fty corresponding to these are not 
necessarily equal, since F<p may be a many-valued function. Thus, if 
(4> — gif = 0, the latent roots of $* are iv'S'i, and the identical equation in 
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Sylvester's Formula. 

§16. Sylvester has given without demonstration the following theorem in 
the Johns Hopkins Univ. Circ. (No. 28, 1884) : 

F<p = XF 9l . (»-fl)(»--ft)- •••(»-&.) . 
(9i — 9%)(9i —9z) (9i — 9u) 

It may be proved as follows : Suppose the latent roots of <£> are all distinct, and 

denote for convenience the left-hand member of the above equation by (2) ; then 

for the a axes of $ we have (I,) p 1 = Fg 1 .p 1 , (X)p z = Fg 2 .p 2 , etc. But by the 

law of congruity and the law of latency, F(<p) p x = Fg 1 .p 1 , F($) p 2 = Fg % .p % , etc., 

for the a axes; consequently (2)pi = F(q>)pi, (X) p 2 = F(<p)p 2 , etc. ; and hence 

for any quantity linear in the axes p = z^ + z 2 p 2 + e ^c., F(<p) p = (2) p. Since 

the latent roots of <£> are all supposed distinct, p may be any quantity whatever 

in the ground of which <£> is a linear unit function. Hence 

i^>=(2). 

As this mode of proof is a verification, if F<p is a many-valued function, 
it is necessary to show that Sylvester's formula gives all possible solutions of 
■4' : =F(<p). Let the latent roots of <p be all distinct; take as latent roots of a 
matrix 4 1 any set of the values of Fg lt Fg z , etc., and as axes corresponding respec- 
tively to them the a linearly independent axes of <£>, p lf p 2 , etc. ; then ^ satisfies the 
equation i^-p = Fty.p for each of these axes, and consequently for any quantity 
in the ground. Evidently m w matrices 4 1 ma y thus be formed if m is the number 
of values of the function Fg, and only the matrices so formed satisfy for all 
values of p the condition 4f = F<p.p, and thus are solutions of ^ = F<P' These 
rrC matrices are, however, all contained in Sylvester's formula. 

By the theory of limits the theorem i5iay be extended to the case where 
the latent roots of $ are not all distinct. 

In the general case in which the latent roots of <?> are all distinct, a simple 
proof of Sylvester's formula may also be derived from the form of F<p when the 
axes uf <£> represent the ground. In §5 it was shown that when the axes of q> 
are linearly independent, whatever the set of quantities chosen to represent the 
ground, there is always a definite matrix w such that 

p = (xg\g 1 fw- 1 ), 

# 3 

g a 
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whence the form of $ being given in terms of any set of quantities, we can very 
simply find the form of F<p in terms of the same set of quantities, for 

F^ = (w\F gi ....0 \w)-\ 
Fg t ....O 



. . . . Fg a 



If <£ is a scalar, any quantity in the ground is an axis, hence the above con- 
siderations showing that Sylvester's formula gives all the solution of ^ = F<p do 
not apply. And on inspection it is evident that the formula does not give the 
non-scalar roots of a scalar. These must therefore be formed by an independent 
investigation. There are evidently two cases to be considered, the roots of the 
matrix zero and the roots of the matrix unity ; from the latter the roots of any 
non-zero matrix may be found. 

Boots of the Matrix Zero. 

§17. It is very evident by §13 that the latent roots of a nilpotent quantity 
are all zero, and, conversely, that a matrix all of whose latent roots are zero is 
nilpotent. Whence it follows that the number of roots of zero with any index 
is infinite, meaning by index of a root of a matrix the least power of the root 
which will reproduce the matrix. Of the square roots of zero in a matrix of 
order a there are o 2 — 1 linearly independent. Thus, if a = 2, the four vids of 
a dual matrix may be expressed linearly in terms of 

(A:A) + (B:B), (A:A) — {B:B)+{A:B) — {B;A), (A:B) and (B:A), 

of which the first is unity and the other three square roots of zero. The propo- 
sition may be shown similarly for any value of o . 

No root of zero can have an index greater than a. For, as Benjamin Peirce 
has shown, there is no linear relation possible between the powers of a nilpotent 
quantity that does not vanish ; but if $ is a root of zero, unless its a m or some 
lower power vanishes, by means of the identical equation, the o th power of q> 
can be expressed linearly in terms of the lower powers, the coefficients not all 
being zero. 

If 4» is an m th root of zero, then by §11 the nullity of successive powers of q> 
increases until the m tb power is reached ; and the nullity of q? is at most twice 
49 
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the nullity of $, the nullity of q> 3 exceeds the nullity of $ 8 by an amount not 
greater than the increment of the nullity of <£> 2 over the nullity of <£> , etc.* Follow- 
ing the reasoning of §9, let the extension of the p linearly independent quantities 
a 1( ag, . . . . a p be the null region of q>; let the extension of the q linearly inde- 
pendent quantities @ x , (3 t , . . . . /3 S be that transformed by $ into the null region 
of q> ; let the extension of the r linearly independent quantities y x , y % , , . . . y r 
be that transformed by <£> into the extension of the /3's, etc.; finally let the s 

linearly independent quantities x x , x % x s be that transformed by 4> into the 

extension annulled by ^> m - z additional to that annulled by lower powers of 4> ; 
and let the extension of the t linearly independent quantities 7^,\, . . . .% t con- 
stitute the remaining extension of the ground. Then these sets of quantities, 
the a's, /2's, etc., together embrace the extension of the ground, and 

<£>«! = 4>«2 =....= <p<x r = , 



$Pi = «n«i + a 21 a 2 +.... + a pl a p , 
$/3 2 = a n a x + a^ag +....+ a p2 a p , 



4>/3 g = a lg a x + a 2s a 3 + 
<?>7i = &uft + &21& + 






etc. ; finally, 



4>7r = &lrfr + &*#■ + + &«r&> 






<^ = \ t K x + h it Xz + 



+ hfi,' 



* Hence certain roots of nilpotent quantities have no representation. Thus let ^ be a matrix of order 
six whose nullity is two, and whose fourth power and no lower power vanishes ; then ^ has no square 
root. For if f 2 = <t> , then, since V has at least simple nullity, the nullity of ty 6 = s is at least six : hence 
S =0. 
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Whence the form now representing <jft is 





a n a x % • • • ■ aiq 

a pl a pZ • • • ■ tt-pg 










b n b n .... b lr 

"21 "22 ••• • b ir 




Og\ o q2 .... b v 






: 






k n & 12 . . . . Jc u 
•hx <hi, • • • • "kt 

K sl fc s2 . . . . Jc st 













In this form all the constituents except those in the rectangles bordered on 
the left and below by heavy lines are zero. The successive heavy lines at 
every other alternate corner intersect with the diagonal, and no one of them 
is longer than the next preceding one. The square at the left-hand upper corner 
and the square . at the lower right-hand corner consist respectively of p* and f 
zeros. 

If <£> is a root of zero, so also is its converse. 
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Boots of the Matrix Unity. 
§18. If 4> is an m th root of unity it is subject to the condition 

*-— i = fo — A0(* — **) .... ($ — a*) = o, 

the Jl's being the scalar m th roots of unity. By §13, if any n < a of the «i Jl's, as 

Jli, X 8 X„ repeated respectively \, \, . . . . k n times (the sum of the &'s 

being o), are chosen as the latent roots of q> ; and if q> — \, $ — a, g , . . . . q> — X w 
have the respective nullities &j, \, . . . . Jc n , then $ satisfies the condition; and 
conversely. It is obvious that the number of the non-scalar roots of unity with any 
index is infinite.* If no lower power of q> than the m th is a scalar, then q> will be 
termed a primitive m th root of unity. The condition necessary ' and sufficient 
that <£> , in addition to being an m th root of unity, shall be a primitive m tlL root, 
is that one at least of its latent roots shall be primitive. 

When the index m<a, and the entire set of the scalar m th roots of unity 
are latent roots of $, then any m successive powers of <|> are linearly inde- 
pendent. For then the identical equation is <p m — 1 = 0; but if there were 
any linear relation between m successive powers of <£, the identical equation 
would be of order rn — 1 , which, by § 2, is impossible if q> has m distinct latent 
roots. 

In the calculus of matrices, as in ordinary algebra, the rn th roots of any 
matrix may in general be obtained by taking any one of its m th roots and mul- 
tiplying it successively by all the to" m th roots of unity of a certain set, namely, 
those that have as axes the axes of the matrix. 

It is obvious that if <£ is an m ih root of unity, so also is its converse, and if q> 
is primitive, so also is $. 

§19. Every matrix whose order is a prime number has u 3 linearly independent, 
primitive, a th roots of unity. For, a being prime, let ty denote the matrix 

*If $ is not a scalar and has all its latent roots distinct, as was stated in 316, i> m (where m is an 
integer) has in general m™ values, obtained by combining each axis of <j> with one of the m th roots of the 
latent root of ^ corresponding to that axis. The non-scalar roots of a scalar may similarly be obtained, 
but since any quantity is the axis of a scalar, these roots' are infinite in number. 
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in which all the constituents are zero except one in each column, and the non- 
zero constituents, forming a broken diagonal, are denoted each by a with sub- 
script indicating the row in which it appears ; the constituent in the X th row 
appears in the first column, that in the (x + l) th row in the second column, and, 
in general, if \jx] denotes the smallest positive residue (modulus a) of as, then the 
constituent in the [x + r — l] th rows appears in the r th column. The constituent 
in the first column I shall term the leading constituent. It is evident, then, if 
(«!, a 2 , .' . . . a„) are the elements of the set upon which q> operates, 

$<Xi = a K a K , 4>o^ =a K + 1 a K + 1 , etc., 
and that the general expression for the <p of any one a r of the a's is 

I. e. <f> applied to the a's advances each by x — 1 places and multiplies it by a 
certain one of the a's. The application of <£> 2 to each of the a's advances it by 
2(x — 1) places, and multiplies it by the product of two of the a's, etc.; finally, 
4> w advances each a by a multiple of a places, i. e., transforms each a into itself 
and multiplies it by the product of all the a's. If, now, X be an imaginary a tb 
root of unity, and the successive a's, a lt a z , . . . . a a are severally put equal to 
1, X, %?, . . . . Jl" -1 , denoting by ft what ^> then becomes, We have 

<o(o.-f 1 ) 

$i («i. <*2» ««.) = *- a (ai, a 2 a,). 

Whence unless a = 2, ft is an a th root of unity ; and since a is prime, it is readily 
seen that no lower power of ft is a scalar. Similarly, denoting by ft what <p becomes 
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when the a's in the same order are put equal to 1 , a, 2 , a, 4 a 2(w-1) , then ft is 

also a primitive o th root of unity ; finally, denoting by ft,.^ what <p becomes 
when for the a's in the same order are substituted the (a — l) th powers of 
1, a,, a 2 , . . . . V~ l , then $>"_! = 1 and also no lower power of ft,-! is a scalar. 
In this way a — 1 matrices are obtained which are primitive a ib roots of unity ; 
besides these, another primitive a th root may be obtained by putting the a's all 
equal to unity, and this matrix may be denoted by ft. These a matrices are 
linearly independent, for otherwise, for some set of values of the x's, other than 
all equal to zero, we should have 

^oft + «i$i + ®$2 +•••• + as<o-ift>-i = , 



Xq + x 1 


+ a» 


+ ••• 


• + a u -i = 0, 


x + a,ccj 


+ tfx 2 


+ ... 


. + ^~ 1 x a _ 1 =0, 


x + a, 2 ^ 


+ %% 


+ ... 


. +a, 2 < w - 1) a; a ,_ 1 = 0, 



which is impossible. Thus any position of the leading constituent other than in 
the first row (i. e. for any value of x other than x = 1), a proper choice of the a's 
gives a linearly independent o th roots of unity ; and thus in this way a {a — 1) 
linearly independent o th roots of unity may be obtained ; for it is evident that 
the matrices obtained for different values of % are linearly independent, being 
expressions in different vids. In addition to these, let a — 1 matrices be formed 
by taking as constituents of the principal diagonal in the first matrix the set of 
scalar 1, a, a, 2 , .... ^" _1 ; in the second, the squares of this set, etc., and finally, 
the (q — l) th powers of this set for the (a — l) th matrix, all the other constituents in 
each matrix being zero. These a — 1 matrices are primitive a th roots of unity, as 
may easily be seen by considering their effect upon the set (a 1( ag, . . . . a M ). 
That they are linearly independent of each other and unity may be shown by 
the above reasoning ; and since they are formed from the vids along the diago- 
nal, which do not occur in the other set of a (»— 1) roots of unity, they can 
have no linear relation with these. Thus, in addition to unity (or rather to a 
primitive scalar u th root of unity), a 8 — 1 other primitive o th root of unity may 
be formed, all of which are linearly independent. 

If a = 2, the four matrices formed in this way are 

(10) (1 0) (01) (0-1 ) 

I 1 I I —1 I 1 1 o I 1 1 oi" 
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Of these the first is unity, the second and third are square roots of unity, and 
the fourth is a square root of — 1 . By multiplying the fourth root by */ — 1 , 
it is of course converted into a square root of unity. Thus a dual matrix may 
also be expressed linearly in terms of four square roots of unity. 

Algebras Analogous to Quaternions. 

§20. It was stated in I, §10, that if i,J, k are any three mutually normal 
unit vectors, any quaternion may be expressed linearly in terms of the four new 

umts 1 +i*f^~l j + kV^l — / + W^l 1 — W^l 

2 2 " 2 2 

which, having the same multiplication table as the four vids of a dual matrix, 

(A:A), (A:B), (B:A), (B:B), 

may, consequently, be regarded as respectively identical with them. This iden- 
tification gives the following values for the ordinary quaternion units, 



1 = (1 o) i = (— V— 1 0) j={ 1) k = ( 0— V — 1) 



|oi| I o V— 1 1 |—i o| |— V— i o| 

The discovery of this form of quaternions, which I have termed the canonical 
form of quaternions (I, §8), as has been stated, is due to Benjamin Peirce ; 
it received its full significance only after the discovery by his son, Charles 
Peirce, of the unlimited' system of quadrates formed from the system of vids 
(A: A), (A:B), etc., when it appeared that quaternions was only the first of this 
system of quadrate algebras, and the identification of quaternions with the 
theory of dual matrices was virtually accomplished. Evidently of all linear 
associative algebras the quadrate algebras form a class which are closely related, 
and consequently are closely analogous to quaternions. In the preceding section 
it was shown that all matrices or quadrates whose order is a prime number may 
be regarded as linear in unity and a z — 1 linearly independent, primitive th 
roots of unity, just as quaternions is an algebra linear in unity and three square 
roots of — 1 , or of unity ; whence the analogy between quaternions and these other 
quadrates extends beyond the quadrate form possessed by each. Moreover, 
these o % — 1 roots of unity are formed from the a 1 vids of the quadrate of order 
0, except for a scalar factor, precisely as the i,j, h of quaternions are formed 
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from the vids of the dual matrix. Thus the quadrate algebra of prime order 
next in order to quaternions (nonions) is that formed from the vids of a triple 
matrix: comparing its eight cube roots of unity with the i,j, k of quaternions, 
we have 



i=[(A'.A)-(B:B)lV-l, j=(A:B)-(B:A), h= [(A:B) + (B-.AftV-l, 

while the nine nonion cube roots are 

(A:A) + h(B:B) + tf(G:G), (A:A) + 7? (B:B) + K(G: G) , 

(A:B) + {B:C) + (C:A) , (A:B) \ + X(B: G) + X (0:A), 

(A:B) + tf{B:G) + ^(C:A) , {A : 0) + {B : A) +(G:B) , 

(A:G) + *{B:A) +X(G:B), (A;G) + 7?{BU) + l(0:B) ;* 

and the representation of unity in the first system is (A: A) + (B:B), and in the 
second (A: A) + (B:B) + (<7: G). It thus appears that there are an infinity of 
algebras exactly analogous to quaternions, namely, those formed from the vids 
of the matrices whose order is prime. I shall now proceed to show that this 
analogy may be still farther extended, and that all the algebras analogous to 
quaternions, and indeed matrices of any order, admit of selective symbols like 
the S and V of quaternions, are resolvable into the product of a versor and a 
tensor, and that there are functions similar to the conjugate of a quaternion, 
equal in number to a — 1 ; if a is the order of the algebra. I shall show this in 
the case of nonions, but the proof will be applicable to any case. 
§21. A quaternion q may, in general, be represented by 

q = Sq + Vq = a + bi, 
where i is a unit vector. In matrix form this is 



q = a(l 0) + 6(V— 1 0) = (a + 6V— 1 0) 



l| I —V— 1| | a — bV— l| 

Hence a + b\^-l = Sq + TVqW^^^l , and a — b\/^~l = Sq— TYq.sf^X 
are the latent roots of q. 

Kq=Sq—Vq=(a — bV=l 0) 

* These units are the converse of those given by the method of the preceding section. The first of 
the series of quadrate algebras analogous to quaternions, nonions, was discovered independently by the 
Peirces and Sylvester. I have in I, ?10, given a short account of this discovery. 
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i. e. the conjugate of q is obtained from q by interchanging its latent roots but 
leaving its axes unchanged. 

Now a nonion, or matrix of the third order, provided its axes are taken to 
represent the ground, may in general be put in the form 



( 


ft 

</ 2 

g 3 


> = ( 


a -\- b -\- c 6 
a + Xb + tfc 



o) 


a + X*b + Xc 


= a \ 


10 

10 


)+M 


1 
A 


)+*( 


10 

o a, 2 o 


), 




1 






o o a 2 




o % 





where Jl is an imaginary cube root, of unity. The first of the matrices in the 
third member is unity, the second is a non-scalar cube root of unity, the third, 
being the square of the second, is also a non-scalar cube root of unity. 
Denoting the second matrix by i, n becomes a-\- bi-\-c$. To select the scalar 
and non-scalar parts of n, the selective symbols S and V may be employed; 
and to discriminate between the first and second parts of the non-scalar portion 
of n the V may be written with subscripts 1 and 2. Employing this notation, 
we have 

n = Sn + V x n + V z n = a-\- bi + ci 2 . 

And since V x n and V 2 n are scalar multiples of non-scalar cube roots of unity, we 
have V(V x n) 3 = V(V z n) s = 0, just as in quaternions Y(Vqf=- 0; and as in qua- 
ternions the last formula gives the important result V(Vq.Vq' + Vq'. Vq) = 0, 
so letting V x n = a 1( V x n' = /? x , V x n" = y x , and V 2 n = a 2 , Y % ri = 9) V 2 n" = y 2 , 
from the two nonion formulae we may obtain in a similar way the two follow- 
ing results : 

^(oi&yi + a i7i/2i + /Vi/i + Piyiai + Yia-fii + 7i/Vi) = , 

and a similar formula in which the a x , (3 X , <y x are replaced by a a , /? 2 , JV 

As in quaternions V 2 q = — I®Vq, so in the algebra df nonions we may 

write b 3 = Vfn = T s V x n , c 3 = PJ» = T s V % n . 

Since the conjugate of q is obtained by interchanging its latent roots, this 
suggests that a cyclic interchange of the latent roots of n, leaving its axes 
50 
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unchanged, should produce a function of n similar to the conjugate of q; and the 
conjugate of n may be defined by the equation 

9s 
9x 

The conjugate of Kn is KKn, which may be written E?n; and evidently 

K*n = ( g, ) = Sn + XV& + KV % n. 

9\ 

92 

These formulae resemble that for the conjugate of q. In quaternions K % = 1 , 
but in nonions we may write K 2 = K~\ and K s = 1 . 

The tensor of a quaternion may be defined as the square root of the product 
of the quaternion and its conjugate. Following the analogy thus suggested, the 
tensor of a nonion n may be defined by the equation 

T 3 n = n.En.K 2 n = a 3 + V +. c 3 — Sale, 
.: T 3 n = S 3 n + T 3 V& + T 3 V % n — ZSn. TV x n. TV % n, 

which is the analogue of !Pq = S*q + 'PYq. It is readily seen that the square 
of the tensor of a quaternion is equal to the product of its latent roots, and thus 
to its content ; and similarly, that the cube of the tensor of a nonion is equal to 
the product of its latent roots, and hence to its content. Whence, since the 
determinant or content of the product of two matrices is equal to the product of 
their contents, the tensor of the product of two nonions is equal to the product 
of their tensors. When a nonion is expressed in terms of unity and eight non- 
scalar cube roots of unity, this proposition gives a theorem analogous to Euler's 
theorem when that is regarded as a theorem relating to the product of two qua- 
ternions ; but the tensor of a nonion is then too complicated an expression to 
give the theorem any interest. 

§22. If i is a. nonion cube root of unity whose latent roots are 1,%,%? 
(yl being an imaginary cube root of unity), and if e denotes the base of the Nape- 
rian logarithms, by Sylvester's formula, 

«« = J w (* + »+i) + te M (tf + a* + *.»)'+ av* (* 2 + a 2 » + a)] . 
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The coefficients of unity, i, and #, in this expression I shall denote by / (0), f x (0) 
and f % (0) ; they are obviously analogous to sin and cos which appear in the 
corresponding expression for e" a in quaternions (where a is a unit vector). This 

g ives «*=/o(0)+ /i(0).*+ /*(0)-* 3 > 

••• £ A9i =/oW+ Vi(0)-* + »."/. W-i 1 , 

«" l =/,(0)+ay 1 (0).» + V.W.*'- 

The second and third expressions are, severally, the first and second conjugates 
of E e \ Since TV* = e W w = 1 , hence 



/oW 8 +/ 1 (0) 3 +/ 8 (0) 8 -3/ o (0)./ 1 (0)./ 2 (0) = l. 

This suggests the corresponding formula cos 2 + sin 2 0=1. The properties of 
the sine and cosine, that cos ( — 0) = cos0, sin ( — 0) = — sin 0, have their 
analogues in these functions ; for from the values for e* w and e Aa< ", it follows that 

MM)=M$), AM = *f 1 (6), AW) = HMO). 

These functions also give rise to a formula in nonions analogous to De 
Moivre's theorem: thus, since 6°V i = e (a+w ', hence 

[/»(«) + /i (*)••* + /.(«) .*][/«(/?) +/i (/?)•» + /» (/?)*] 

= /o(« + /?)+/i (« + £).*+/>(a + /?).** 

Whence arise formulae for the functions of the sum of two arguments analogous 
to the formulae for the cosine and sine of the sum of two angles, namely, 

/.(a + 0)=/.(a)./oGS) +/i («)/,(£) +/.(a)/iG3), 
/i(« + /?) =M*).fi(P) +fi(*).fo(P) +/.(«)/. 03), 
/■ (« + /»)= /o («) •/> 0?) + /■ (a) ./„ 08) + /i (a) .A 08) • 

It is now easy to perceive that, — just as for any quaternion q, for a proper 
value of 0, 

q = e>°«*r+t«->g*. tfr 9 _ ^ ( cog q + gin 0> jjy^ . 

so for any nonion n for proper values of and yj , 

n = e losFn + ii + ^= Tn(f 6 +f i e.i+f i e^)(f y l +/ lVI ^ +/ i y,.i) 

= Tn(foW + l)+fiV + v)i+M9 + *)*)- 
We have therefore 

n = Tn. Un, 

n.ri= Tn.Tn'.Un.Un'. 
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§23. Another analogy exists between quaternions and the quadrates of 

prime order, namely, that just as quaternions is linear in unity and two square 

roots of unity and their product, so the matrices of order a constitute an algebra 

linear in unity and two o th roots of unity, their powers and products. Thus in 

the case of nonions, as Sylvester has shown,* if 7. is an imaginary cube root 

of unity, and 

* = ( 1) j — ( 1) 





7? 



7? 






a, 



then ij — 7?ji and ji = Tiij, while the products formed from the two sets (1 , i, t 2 ) 
and (I,/,/ 2 ) give all the nine units of the octanomial form of nonions. More- 
over, I find that just as there are an infinite number of systems of three mutually 
normal unit vectors, so there is an infinite number of systems of «ight cube roots 
of unity similar to the system formed from i and j and their products. The 
general expressions for these two new cube roots i x and j\, from which the new 
system is to be formed, are 



(- 









1 


A, 
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„— 1 



) *=(. 
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7? 
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- 1 ) 



where \m j = ; If the axes of i are (p x , p 2 , p 3 ) and the axes of/ are (<s x , <r 2 , <r 3 ), 
then the axes of i x and j\ are respectively (vtp lf wp z , wp 3 ) and (wa lt wa % , md z ). 
The condition that from i x and j\ an octanomial system shall be formed similar 
to that formed from i and/, is that i x and j\ shall have as latent roots the three 
cube roots of unity, and that the axes of ix and j\ shall be related in the same 
way as those of i and/. If (a lt <x 2 , a 3 ) represent the ground 



p x = ai + Aa 2 + a 3 , 
P2 = «i + <*2 + Xa 3 , 
p 3 = Xa] + a 2 + a 3 , 



Oi = % + tfa 2 + a 3 , 
2 = tfdi + otg + a 3 , 

$3 = O-l + «2 + ^ 2 «3 • 



With respect to the general case, let i be a primitive co th root of unity 
formed according to the method of §19, whose leading constituent is in the k th 
place, and the coefficients of whose successive vids (beginning with that in the 
first row) are 1, 7,, 7?, . . . . 7J"~ 1 (% being an imaginary cj th root of unity); and 



* Johns Hopkins University Circulars, No. 17 (Aug. 1882). 
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let/ be formed from the same vids with the coefficients oN replaced respectively 
by their squares ; then ij — %"-* +1 ji and ji=.,h k - l ij. If k = 2, ij = V~ l ji and 
ji = Mj. 

Quadrate Algebras whose Order is not a Prime Number. 

§24. With regard to the quadrates whose order is not a prime number, it 
may be shown also by a method similar to the method of §19, that they possess 
o 2 — 1 linearly independent non-scalar co th roots of unity, formed from their vids 
in the same manner in which the i,j, k of quaternions are formed from the vids 
of a matrix of order two. However, of these 6> 2 — 1 a tb roots, only (a + l).v (a) 
(where t (a) denotes the totient of «.) are primitive, namely, t (&>) o th roots formed 
from the vids along the diagonal, and a.? (a) others formed from the non-sym- 
metric vids. The latter consist of those roots in which the leading constituent 
is in a row whose order less one is prime to a. The roots in which the leading 
constituent is in a row whose order less one has with « the greatest common 
divisor 8, other than unity, are either (o : 8) th roots of unity, or such to a factor 
pres. Of the th roots of unity formed from the vids along the diagonal by taking 
as their coefficients powers of the scalar th roots of unity, 1, ?l, 2?, . . . . Jl w-1 , 
only those are primitive 6> th roots of which the power taken of the set of 3,'s is 
an integer prime to a. 

However, all these quadrates are analogous to quaternions in admitting of 
selective symbols, in having functions analogous to the conjugate in quaternions, 
and in that any expression in them is resolvable into the product of a tensor 
and a versor. This may be proved in precisely the same way in which in §21 
and §22 it is shown to be true for the quadrate algebra of order three. 

The quadrate algebras whose order is not a prime number are compounds 
of other algebras, by which I mean that they are linear in the products of the 
units of these algebras : if a = «'. a", then the quadrate algebra of order a is 
linear in the products of the units of the quadrate algebras of order a' and a", 
the units of each of these systems being regarded as commutative with those of 
the other. Thus, the algebra formed from the vids of a matrix of order four is 
linear in the products of two quaternion sets, the i,j, k of the one set being 
commutative with the i, j, k of the other set ; the algebra formed from the vids 
of a matrix of order six is a compound of a quaternion set and a nonion set, 
the units of the quaternion set being commutative with those of the nonion 
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set. It is obviously sufficient to prove this for any form of the several algebras, 
and for this purpose I choose the canonical form, and I shall show that the vids of 
any matrix of order a' combined with those of any matrix of order q" (the vids 
of the first set being regarded as commutative with those of the second set), give 
an algebra whose multiplication table is the same as that formed from the vids 
of a matrix of order a = o'.o". I shall, however, first illustrate this theorem by 
the case of the matrix of order four. Considering the two quadrate systems 

(A.-.A,) (A.-.A,) (B.-.B,) (B X ;B % ), 
(At-.AJ (A 2 :A 2 ) (B.-.B,) {B 2 -.B % ), 

which will be regarded as commutative, evidently the complete system formed 
by their products will consist of sixteen linearly independent units which may 
be arranged as follows : 



(A,: A^B.-.B,) (A.-.A^B.-.B,) 
(A.-.A^B.-.B,) {A,:A,){B X :B X ) 



{A l :A 1 ){B t :B0 (A 1 :A i )(B,:B 1 ) 
(A i :A 1 )(B,:B 1 ) (A 2 : A 2 )(B Z : B,) 



{A 1 :A 1 )(B l :B t ) (A.-.A^B.-.B,) 
(A.U^B.-.B,) (Az-.A^-.Bz) 



{A.-.A^B.-.B,) (A.-.A^B.-.B,) 
{A.-.A^B.-.B,) (A % :A,)(B % :B 2 ) 



Those compound vids in the upper left-hand group consist of the first set, each 
multiplied by (B^.B^); those in the other groups also of the first set multiplied 
respectively by (B X :B^, {B Z :B^) and (B 2 :B 2 ). According to this scheme the 
product of the vid (A x :A y ) from the first system and the vid (B r :B ( ) in the 
second system occupies in the resulting system a position which may be denoted 
by (x + a'r — 1 , y + dt — 1) , where, as before, the first number denotes the 
row and the second the column; and in this case a' = a" = 2. On trial it will 
be found that this compound system has the same multiplication table as that of 
the system of vids of a matrix of order four, the vids in the compound system 
corresponding to those of the quadrate system of order four which occupy the 
same place. It should be observed that if I and i 7 are vids of the first quadrate 
system and J and J' vids of the second, that the product of the compound vids 
I.J and I'. J' is zero, unless I.I' =f= and J.J' zf= o ; since the I's and J's are 
commutative. The proof in the general case may be accomplished as follows : 
Following the same scheme of arrangement of the compound vids formed from 
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the product of the vids of two quadrate systems of order ©' and a" respectively, 
we have 



IJ= (A x :A y )(B r :B s ) =(x+ o 'r — l , y + at e - 1 ) = K, 
• I f J f =(A y ,:A s )(B 8> :B t ) = (y> + g>V=1, z + o'V^l) = K>. 

The product of IJ and I' J' is zero, unless y' =■ y and s 1 = s, when we have 

7/.J'«7'=(^:A)(^:^), 



which in the compound system will occupy the (x + a'r — l) th row and 
(a -+- alt — l) th column, and thus is represented by (x + a'r — 1, z-\-'a/t — 1). 
The necessary and sufficient condition, however, that the compound system shall 
have the same law of multiplication as that of the algebra formed from the vids 
of a matrix of order a = «'.«" is that KK' shall be zero unless 



y' + id's' — 1 = y + a's — 1 , 
in which case we must have 



KK' = (x + o'r — l, z + a't — 1). 



But since y' < «', y < a', if y 1 + «'«' — 1 =z y + «'« — 1 , then y' = y and s' = s. 
Hence the condition is fulfilled. 

§25. From the last section it follows that the matrix of order a = 2 m is a 
compound of m quaternion algebras which do not interfere, i. e. the units of each 
set are commutative with those of the other sets. I shall term such an algebra 
the tn-way quaternion algebra. The system of m-way quaternion algebras has 
already been considered by Clifford (this Journal, Vol. I). Clifford, however, 
approached the subject from an entirely different point of view. He starts with 
n "elementary units" i x i % . . . i n whose multiplication with each other is' polar, and 
which are such that the square of each is — 1 . The 2" products, of order to n , of 
these elementary units, are linearly independent ; and Clifford shows that the 
products of even order, 2 M_1 in number (which may be obtained from all combi- 
nations of the binary products ifa, ifo, etc), form an algebra or system by them- 
selves, which he terms the n-way algebra. If n =■ 2m + 1 , Clifford further 
proves that " the n-way algebra is a compound of m sets of quaternions which do 
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not interfere ; " and if n = 2m , he shows that the n-w&y algebra is a compound 
of tn quaternion sets which do not interfere, and the algebra (1 , g) , where e 2 = 1 
and g is commutative with each of the quaternion sets. At present we are only 
concerned with the n-way algebra when n is odd ; it consists, we have seen, of 
2" -1 = 2 2m linearly independent units, and is a compound of rn quaternion sets 
which do not interfere. Hence Clifford's n (,= 2m + l)-way algebra is that 
formed from the vids of a matrix of order 2 m , or is the m-way quaternion algebra. 
Commenting on the surprising fact that sets of quaternions should appear as 
the simplest form of an algebra which at first sight is so far from suggesting 
Hamilton's system, Clifford says that "thus it appears that quaternions is the 
last word of algebra to geometry." It is still more surprising that quaternions 
should be so prominent in the theory of matrices, and in a sense embrace the 
whole subject; for since the theory of matrices of any order is included in the 
theory of matrices of higher order, and since however great a number may be, a 
power of two may be obtained which is still greater, it follows that the theory 
of matrices of any order is included in the theory of the combination of a certain 
finite number of quaternion sets which do not interfere. 

With regard to Clifford's geometrical algebras I am able to show that the 
entire system formed from the combinations of all orders of the n elementally 
units, which I term the w-fold algebra, is identical with the (n — l)-way algebra ; 
and when n= 2m, by an argument similar to that employed in the case of 
bi-quaternions by Benjamin Peirce, that the n-way algebra is a compound of two 
groups of m sets of quaternions which do not interfere, such that every set of 
either group is nil-facient and nil-faciend with any set of the second group. 

It is clear that any quantity in Clifford's n-fold algebra may be expressed as 
a sum of products of expressions a = Xa x i v ft = 2&i*i, etc., linear in the elemen- 
tary units; the a, ft, etc., are of course such that a 2 = — 2a*, ft 2 = — %b\; 
whence, on account of the obvious analogy, they may be termed vectors. In 
virtue of the relation stated above which a matrix of any order holds to those 
of higher order, and the consequent inclusion of the theory of matrices of any 
order in that of sets of commutative quaternions, the theorem follows that the 
theory of matrices is the theory of expressions which are sums of products of 
quantities a, ft, etc., whose squares are negative scalars. 
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The Identical Equation. 

§ 26. The representation of any matrix as linear in unity and the o — 1 
power of a non-scalar a tb root of unity gives a very simple demonstration of the 
identical equation in the case of quaternions and of nonions. Thus, 

(q—Sqf= V*q = -T*Vq, 
.'. <f—2Sq.q+ T*q=0, 

(n — Snf = ( Vfl + V % n)' = {hi + d*) 3 
= b 3 + c 3 + Sbc (bi + e?) 
= T 3 Vji. + T 3 V,n + 3 TV x n TV 2 n (n — Sn) , 

.: n 3 — 2,Sn.n % + 3 (S*n + TV,n. TV z n.Sq) n— T 3 n=0. 

By the following method the identical equation is proved immediately for 
any case. It, however, depends upon the lemma that any symmetric function of 
a conjugate matrices is a scalar, and is the same symmetric function of their latent 
roots. The lemma is very readily proved by considering the matrix in the form 

9i( a i :a i) + 9% (oa^oa) + •••• + &. (««:«»), 

when the conjugates appear as 

ft(oi:oi) + 9>(a*'-<h) + • ■ ■ ■ + 9\ K : ««,) . 

etc. Returning to the identical equation, we have for a = 3 , 

= (n — n)(n — Kn)(n — K % n) 
= n 3 — (n + Kn + IPn)n 2 + {Kn.K^n + JPn.n + n.Kn)n — n.Kn.IPn 
= n 3 — (ft + ft + ft) n* + (ftft + ftft + ftft) n — ftftft . 

The same proof applies in any case. 



51. 
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POSTSCRIPT. 

When the above paper was written I did not know of Mr. Bucheim's 
paper on the Theory of Matrices in the Proc. London Math. Soc. (Vol. 16). Mr. 
Bucheim's proofs of Sylvester's law of nullity, of the identical equation, and (by 
implication) of the law of latency and what I have termed the law of congruity 
of the axes, are substantially the same as mine. I have further completed the 
investigation of the corollary of the law of nullity, have thoroughly treated the 
roots of unity and of zero, and have shown that there is an infinity of linear 
algebras, constituting a sequel to quaternions and nonions, whose laws of combi- 
nation are equivalent to those of matrices of prime order. I have also shown 
that the laws of combination of matrices of composite order are identical with 
those of algebras whose units are the products of the units of linear algebras 
analogous to quaternions and nonions. I have treated the whole subject more 
in detail and more systematically than Mr. Bucheim. 



